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Abstract 

We study the imprint of new particles on the primordial cosmological fluctuations. 
New particles with masses comparable to the Hubble scale produce a distinctive 
signature on the non-gaussianities. This feature arises in the squeezed limit of the 
correlation functions of primordial fluctuations. It consists of particular power law, 
or oscillatory, behavior that contains information about the masses of new particles. 
There is an angular dependence that gives information about the spin. We also 
have a relative phase that crucially depends on the quantum mechanical nature of 
the fluctuations and can be viewed as arising from the interference between two 
processes. While some of these features were noted before in the context of specific 
inflationary scenarios, here we give a general description emphasizing the role of 
symmetries in determining the final result. 



1 Introduction 


Inflation seems to be the highest energy observable natural process. The Hubble scale, H, 
during inflation could be as high as 10^^ GeV. This is much larger than any energy we 
can expect to achieve with particle accelerators in the foreseeable future. It is therefore 
interesting to understand how to extract information about the laws of physics at those 
scales. We can view inflation as a particle accelerator. The physics during inflation is very 
weakly coupled, as compared to usual particle physics. On the other hand going between 
the physics during inflation and astronomical observations we can do today is somewhat 
involved. It is similar to the process of going between the underlying high energy process 
and the detector signals in a hadron collider, except that in cosmology the hnal state 
involves galaxies instead of hadrons. In a hadron collider we look at jets or patterns of 
energy deposition on the detector. In cosmology we similarly look for patterns in the 
distribution of galaxies or in the cosmic microwave background. These patterns appear to 
arise from a single fluctuation mode, the curvature or adiabatic mode conventionally called 
C. We are going to take for granted the very delicate process of going from present day 
astrophysical observations and the primordial shape of (. In order words, we are going to 
imagine that astrophysicists tell us what ({x) is at the end of inflation, before it reenters 
the horizon. What can we learn from this?. How do we extract the physics during inflation 
from this? 

The discipline of collider physics involves going from the direct collider observables 
to the underlying lagrangian of the theory. One of the simplest questions one can ask is 
how to recognize the presence of new particles. In colliders the answer to this question 
is simple, one collects groups of particles (pairs, for example) and one plots the invariant 
mass. If a bump is seen in this distribution, one says that there is a new particle. One 
can also look at the angular distribution of the particles and read off the spin of the new 
particle. 

In this paper we answer a similar question for the “cosmological collider”. Namely, we 
ask: how do we recognize the presence of new particles during inflation ? The simplest 
version of this question is whether we had any new light or massless helds during inflation. 
These are models which have more than one held and can have rather striking signatures, 
such as isocurvature huctuations or violations of the inhationary consistency condition in 
the squeezed limit of the three point function. There is an extensive literature on this 
case, see DEI for reviews. We will concentrate on the case where we have particles whose 
mass is comparable to the Hubble scale. Such particles can be produced by quantum 
huctuations during inhation and they can then decay to ordinary inhatons. This process 
leaves a statistical imprint on the spectrum of primordial huctuations. More precisely, this 
process gives rise to non-gaussianities, see [3] for a review. A particular model of this kind, 
“quasi-single-held inhation”, was studied in mum and some of the features we discuss 
below were described in jS]. 

Of course, even self interactions of the inhaton, or interactions of the inhaton with 
gravity, can give rise to non-gaussianities. We want to know: 
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Which specific non-gaussianities are a signature of new particles during inflation, 
as opposed to signatures that arise due to inflaton self interactions. 

Now, if we think about very massive particles, m 3> hf, then we can integrate them 
out and they produce new terms in the effective lagrangian for the light helds. Since we do 
not know the original lagrangian, it is clear that we are not going to discover them in this 
way. For masses of order H, m ^ H, the situation is different, because we can produce 
the particles giving rise to non-local effects which cannot be mimicked by changing the 
interaction lagrangian of the inflaton. The fact that this is an interesting question was 
addressed in |l El El El Eli- 



Figure 1: (a) The three momenta of the three point function form a closed triangle. In 

the squeezed limit, one of the sides, fciong, is much smaller than the other two, fcshort- (b) 
In this limit we can also dehne the angle 7 between the short and long mode momentum 
vectors, (c) In position space we have two insertions at a short distance from each other, 
associated to /cshort and one at a longer distance, fciong- We are interested in considering 
the effects of massive helds, a, that can decay into pairs of inhatons. In an inhationary 
background we can replace an inhaton huctuation by the classical fio background, so that 
we get a contribution to the three point function. 


The large masses, this ehect is suppressed by which is why it isn’t captured 

by the ehective held theory which is an expansion in powers of {H/m). 

The simplest non-gaussian observable is the three point function Fourier 

space. Translation invariance implies that the momenta form a closed triangle. The near 
scale invariance of the huctuations implies that this is a function of the shape of the 
triangle. In particular we can form the ratio between the smallest side and the largest side 
^long/^short- See hgurej^a). As a function of this ratio the three point function can display 
interesting power law behavior for small values of this ratio. This is called the “squeezed” 
limit. We hnd 
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where Aj can be real or complex, Wi are some coefficients and e is a slow roll parameter 
(see m- This form of the correlator is a consequence of the slightly broken conformal 
symmetry of the late time wavefunction of the inflationary universe. 
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Such power law behavior is characteristic of conformal held theories, where the Aj 
that occur are the (anomalous) dimensions of operators in the theory. Here we have the 
same situation because the de-sitter isometries act on late time expectation values in the 
same way as the conformal group in one less dimension. In inhation, these approximate 
symmetries govern the behavior in the squeezed limit. 


The actual values of Aj that occur in (1.1) give us the desired information about the 
spectrum of new particles. Namely, a single new particle of mass m > ^ and spin s gives 
rise to a pair of terms with [S] 
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Ps{cos9) (1.3) 


A single particle gives rise to two dimensions, A±, which are complex conjugates of each 
other, which is good, so that we can get a real answer for the correlatoiQ An interesting 
fact is that there is an interesting calculable phase in the coefficients w± = |rc±|e- 


± 2(5 


that appear in (1.1). This phase, 5(/i), depends only on the mass of the particle. The 
oscillatory behavior as a function of the ratio of scales is a quantum effect. It arises from 
the quantum interference between two processes. One process is the ordinary gaussian 
quantum evolution of the inflaton. The other is the creation of a pair of massive particles 
that subsequently decay to inflaton modes. We know that we can assign a characteristic 
time during inflation to a given momentum k. It is the time when this mode crossed the 
horizon during inflation. The mode /ciong crossed the horizon earlier in than the mode fcshort- 
The time measured in number of e-folds between these two events is N = log(fcshort/fciong)- 
What we see in (1.3) is the behavior of the amplitude of the wavefunction of the massive 
particle as a function of time, or N. The wavefunction oscillates as we expect for a 


massive particle. In addition the 3/2 in (1.2) gives rise to the expected dilution factor 
due to the expansion of the universe. Namely the square of the wavefunction should go 
like ~ 1/Volume. The fact that we have effects involving the amplitude of the 

wavefunction of the massive particle as opposed to its square is due to the fact that this is 
a quantum interference effect. Another related fact is the dependence of the overall factor 


in (1.3) on the mass. Keeping the coupling between the inflaton and the helds hxed, this 


factor goes as This should be compared to the thermal factor going like This 

last factor is proportional to the probability for creating a pair of massive particles. The 
fact that we get e~'^^ is again due to the fact that we are seeing an interference effect, 
sensitive to the amplitude and containing phase information. 


^ The correlator is real for the following reason, (j(^) is real in position space. Therefore, in momentum 
space it obeys C| = C_fe- H we complex conjugate a correlator, we get the same answer as reversing the 
sign of all momenta. For a two or three point function, we can relate the correlator with reversed mometa 
to the original one by a two dimensional rotation, so that they should have the same value. 
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Needless to say, the observation of non-gaussianities with such oscillations and with 
the precisely predicted phase would be a striking signature of the quantum mechanical 
nature of cosmological fluctuations. We can view this as a kind of “cosmological double 
slit experiment”. The precise prospects for actually measuring this will not be evaluated 
precisely in this paper. However, let us make some comments. For natural values of the 
couplings this is harder to measure than the standard inflationary three point function 
[To] , which does not appear measurable in the near future. However, it is possible to 
have larger couplings, so that it would be interesting to look for such templates in the 
data. Of course, detection of non-gaussianity is limited, in an irreducible fashion, by the 
cosmic variance noise (the fact that we observe only one universe). The squeezed region 
we are studying has a smaller signal to noise ratio than that of a more equilateral triangle. 
Therefore it is likely that one would hrst hud evidence for non-gaussianity and then, by a 
more sensitive measurement, hnd the detailed behavior in the squeezed limit. Our reason 
for emphasizing the squeezed limit is not that this is the easiest to measure, the reason 
is that this limit contains very direct information about the spectrum of particles during 
inflation. This is somewhat similar to the fact that at the LHC we concentrate on events 
with large transverse momentum. This is small subset of the events, but is the subset 
that can display the existence of new particles a direct way. Note that the standard non- 
gaussianity for the single held model HD] can be understood as arising from the effects 
of graviton exchange, as it was explicitly shown in HU. This also shows that the effects 
we are after are expected to be smaller than the ones of the standard non-gaussianity, 
at least if the couplings are all suppresed by the Planck scale. In other words, we can 
consider couplings to the inhaton of the form A f(V(/>)^a, where a is the massive held. 
We generically expect such couplings to be present with A oc l/Mpi. But of course, the 
couplings could be larger. Special models with possibly larger values of non-gaussianity 
due to helds of masses of order Hubble were discussed in H El m E] under the name 
of “quasi-single-held inhation”. In those papers the relevance of the squeezed limit for 
reading oh the presence of new particles was emphasized. 

If the particle is light, with m < 3iL/2, we still get two powers as in (1.2) (1.3) except 
that now A± are both real and A_, being smaller, dominates the squeezed limit. Here 
also the relative coefficient between the two terms is completely hxed by the value of the 
mass. 

It is also possible to study the spin of the particle. If the massive particle has spin, s, 
then we get a factor of a Legendre polynomial, Ps{cos6), in front of (1.3), where cosy is 
the angle between the small momentum fciong and the large one fcshort- 

So far we have discussed the ehects of single particles which can directly couple to 
the inflaton via J (V0)^ct. There are situations where such a coupling is not allowed by a 
gauge symmetry. For example we can have a coupling of the form A where h 

is the Higgs held. If the electroweak symmetry is unbroken, we can only couple to a pair 
of Higgs helds, not to a single Higgs held. See hgure|^d). This situation also leads to an 
expansion of the form (1.1), except that now the leading dimensions that appear rehect 
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the presence of a pair of particles and have the form 


A = 3 + 2i/i , A = 3 — 2i/i , or A = 3 


(1.4) 


where /i is given in terms of the mass of the Higgs held by (1.2). This gives rise to more 
rapidly decaying terms in the squeezed limit, with a dilution factor of 3 instead of 3/2. In 
a similar way one can consider couplings to the Yang Mills helds or pairs of fermions. 

A prevailing attitude in the study of inhationary models is the idea that we have a single 
light held. However, we know that now we do not have a single light held. Where were all 
the helds we have now during inhation?. They could be very massive and unobservable, 
or some of them could be massless (the gauge helds or the fermions) and have observable 
consequences in the squeezed limit. In this case, the natural values of the coupling are 
smaller, ^ f (V(/>)^h'^h, but they could also be bigger!. 


More generally, the powers appearing in (1.1) seem to correspond to all the quasinormal 


modes in the static patch of de-Sitter space. Therefore, in principle, the squeezed limit 
contains the full information about the spectrum of the theory. Of course, in practice we 
are limited by cosmic variance. 





(a) 


(b) 



(c) 



Figure 2: (a) Two point function of massive helds (b) Three point functions of two 

massless helds and a massive held, (c) Four point function of massless helds with an 
exchange of a massive held. We can view this diagram as the amplitude to create a pair of 
massive helds which then decay into inhatons. (d) Contributions from Higgs held exchange, 
or the exchange of any other gauge non-invariant held. 


This paper is organized as follows. We start with a discussion of the two point functions 
of massive helds. See hgurej^a). We then discuss three point functions of two massless 
helds and a massive held, see hgurej^b). We then discuss the ehects of a massive held on 
a massless held in de-Sitter, see hgurej^c). We study the de-Sitter case hrst because the 
result will be largely determined by the de-Sitter symmetries. These de-Sitter symmetries 
act like conformal symmetries on the late time three dimensional slices of the universe. 
In de-Sitter, the main ehect of the massive held appears in the four point function of 
the massless held. In this case, the interesting limit is analogous to an operator product 
expansion were we consider the four point correlator in a kinematic conhguration where, in 
position space, the four points form two pairs which are far from each other, see hgurej^c). 
In this limit the answer will be related to products of three point functions. It turns out 
that the formulas are simplest when we consider a conformally coupled held, = 2H^, 
instead of a massless held. Therefore we analyze this case hrst, and then the massless case 
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is simply obtained from the conformally coupled case. We also discuss general aspects of 
the analytic structure of correlators in momentum space. It is interesting to understand 
where the correlator is smooth and where it is singular. As a toy problem we also hrst 
briefly discuss equal time flat space correlators. 

We then go to the inflationary case, which is related to the case we discussed above 
by replacing one of the external massless helds by the classical inflaton solution 4>o{t), see 
figure [^c). 

We treat both the case of an intermediate massive scalar held and also a massive held 
with spin. The discussion of the spinning case is a bit more technical and the reader 
might want to skip it in a hrst reading. There are two physical motivations for including 
helds with spin. There may be extra dimensions with size not far from the Hubble scale 
during inhation, in which case we can expect the inhation to couple to graviton Kaluza- 
Klein modes. It is also conceivable to have a very “stringy” model of inhation where the 
string scale is comparable to the Hubble scale, with a tiny string coupling. We do not 
have an explicit model of this type, but that could be due simply to a “looking under the 
lamppost’” ehect, namely an explicit construction of such a model lies somewhat outside 
the range of our technical tools (for some attempts in this direction see [I2])- Of course, 
if we were to detect the contribution of a spin 4 state in the non-gaussianity, it would be 
a strong indication of string theory during inhation, since we suspect that a structure like 
that of string theory follows when we have weakly interacting particles with spin s > 2. 


2 Generalities 


In the hrst part of this paper we consider quantum helds propagating on a hxed de-Sitter 
space. In other words, we do not have dynamical gravity. We write the metric as 


ds'^ = 


—dill‘s + dx'^ 


(2.6) 


and we set the Hubble scale to one, H = 1. We will restore it later. 

De Sitter has a large group of isometries. Some are manifest in (2.5). The two that 
are not so manifest correspond to the Killing vectors 


D ; -r]drj- 


b.K : — 2{h.x)7jdn — 2(h.x)x.ds — {rf — \x\^)h.d^ (2.6) 


If we have a quantum held theory in de Sitter, in a de Sitter invariant vacuum, then the 
correlation functions are invariant under the action of these isometries. 

We will be often interested in the late time limit of correlators where some points 
approach the asymptotic future at p = 0 keeping x hxed. In many circumstances one can 
argue that the late time dependence of the operator is a sum of simple exponentials of 
proper time of the form or This is related to the quasi normal mode behavior in 
the static patch. In other words, inside a correlation function, a four dimensional operator 
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can be expanded as 


(2.7) 


{A} 

where we snni over varions valnes of A. Here 0{x) behave as three dimensional operators. 
If (f) is an elementary free field we snm over two valnes of A, as we will discnss later. 

Under these circnmstances, we hnd that the above Killing vectors lead to an action of 
the form 


D : — A — x'^d^i , b.K : — 2A(b.x) — 2{h.x)x.ds + \x\%.ds (2.8) 


on the operators Oa(x). Here we have nsed that r/ —)■ 0. This is precisely the form of the 
generator of the conformal gronp in three dimensions for an operator of dimension A. In 
conclnsion, if the late time expectation valnes of valnes of operators in de-Sitter contain 
simple exponentials of proper time, then the coefficients of these exponentials behave as 
primary operators of a conformal held theory with dimension A. This is a conseqnence of 
de Sitter dynamics and we are not nsing any hypothetical dnality her^ Here 0/^{x) are 
standard operators acting on the Hilbert space of the qnantnm held theory that lives on 
de Sitter space. 

The fact that de-Sitter gives a scale invariant spectrnm of hnctnations is well known. It 
is also trne that the spectrnm is conformal invariant. The extra symmetries are the special 
conformal generators Kj dehned in (2.8). Of conrse the manifest translation symmetries 
x* —)■ X* -|- e* imply that it is convenient to go to Fonrier space. Then the conseqnence of 
the symmetry is momentnm conservation ^ /cj = 0. Note that these are three dimensional 
vectors. Any correlator in fonrier space has the form 


{% ■ • • 0£.> = fe)(OE, ... 0-J 


(2.9) 


where the prime indicates that we have extracted the delta fnnction prefactor associated to 
momentnm conservation. When we compnte scattering amplitndes in hat space we have a 
similar delta fnnction except that we also have an energy conservation delta fnnction. De 
Sitter correlators have an interesting analytic strnctnre which we will discnss in this paper. 
One featnre that we wonld like to point ont here is that pertnrbative correlators have a 
singnlarity when kt —)■ 0, where = |/ci| -|- |/c 2 | -l- ■ ■ ■ |A:„|. Of conrse this cannot happen 
for non zero real valnes of the fcj. However it can happen if we analytically continne the 
ki. The coefficient of this singnlarity is related to the high energy limit of the hat space 
scattering amplitnde 

^ + C.C. , for ^ 0 (2.10) 


^Of course, these properties are necessary for a dS/CFT duality |T31[T3] and were discussed previously 
in that context. In AdS we have a similar property, except that we are computing a special transition 
amplitude where, for elementary fields, the boundary conditions imply that one of the operators vanishes. 
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Here A! is the high energy limit of the amplitude, with all four energy momentum conserva¬ 
tion 6 functions stripped off. The amplitude arises when we evaluate the all -|- contribution 
of the in-in (or Keldysh) contour and the c.c. part comes from the all — part of the contour. 
The singularity at fct ~ 0 comes from a region where all interactions of the perturbative 
diagram are happening at early times, as in hgure (See also P!). 



Figure 3: The singularity at —)■ 0 arises when the whole diagram occurs at very early 
times. This translates into very short distances where we can make the flat space and high 
energy approximation, obtaining (2.10). 


When we perform a cosmological measurement we are limited by cosmic variance. This 
can be alleviated by doing the same measurement in different patches and averaging over 
the patches. When we compute the correlators in Fourier space we are typically performing 
this spatial average already. This is one reason that it is convenient to think in terms of 
Fourier space. If we had the exact conformal symmetry, it is sometimes easier to go to 
Mellin space (see e.g. CS]), and it would be interesting to explore it in the de-Sitter 
context. Here we will stick with fourier space because we will also use it when we go 
away from the exact conformal symmetry in the context of inflation. It is convenient to 
know how the generators of dilatations and conformal transformations act in Fourier space. 


Doing a simple Fourier transformation of the action of (2.8) we hnd 


D: -(A-3) + fc.afc, (2.11) 

b.K : - (A - 3)2b.dk - [{b.k){dk.dk) - 2{k.dk)(b.dk)] (2.12) 

These are the generators acting on each scalar operatoiQ The object that annihilates 
the amplitude {O ■ ■ ■ O)' is the sum of each generator acting on each operator. For the 
dila tation we have that —3 + 'Y^Di annihillates the stripped correlator, where Di is given 
by (2.11) with k ^ ki and A —)■ Aj where i = I,-- - ,n runs over the n operators in 
the correlator. For special conformal transformations 'Yhi b-Ki annihilates the stripped 
correlator. 


3 Two point functions in de Sitter 


The two point function for a scalar held in de-Sitter space depends only on the proper 
distance between the two points. It can be obtained by imposing the equation of motion 


See (A.1571 for the case with spin. 











on the time ordered expectation value 






X) 


=5‘^(x — x') 


(3.13) 


In addition we need to impose that the function is regular at the antipodal point. This 
completely determines the solution 

(0(a:)0(a;')) = ^^2 + + i/i, ^ - ip, 2; 1 --) (3.14) 

(^47rj^ 2 2 w 


w = A 


p = 


r]r] 


— {rj — ri'Y + IT — a;'|2 + ie 

m? 9 

1P~4 


( 3 . 15 ) 

( 3 . 16 ) 


The antipodal point corresponds to tc = 1 and we see that the function is regular there. 


The coincident point singularity is at tc = cx) and the prefactor in (3.14) ensures that 


(3.13) is obeyed. For heavy enough helds p is real. However, these formulas are also valid 


for purely imaginary p, which is the case for light helds. We will be often interested in the 
limit when we keep the comoving distance |x — T'| hxed and we send p, p' —?■ 0. In this 


limit. It is easy to see from (3.14) that the correlator behaves as 


(0(p,T)0(p',O)) 


47J-5/2 


pp 

|T|2 


l+l/J 


r(-^h)r(- Tip) + 


pp 

|Th 


f-iAi 


r(^/i)r(- -fp) 


+■ ■ ■ 


(3.17) 

For real p these are the leading terms. We see that we get two simple exponentials of proper 
time with p = —e“*. As we discussed above, each coefficient of the exponential is 

expected to behave as correlators of operators with the corresponding scaling dimensions. 
This is indeed the case in (3.17). Note that each held gives rise to two operators with 
dimensions 

3 ~ 3 

A = - + zp , and A = 3 — A = - — zp (3.18) 

These operators simply correspond to the annihilation and creation operators with re¬ 
spect to the out vacuum, respectively. The factor p^ ~ is the expected time 


dependence multiplying the out vacuum annihilation operator. We can contrast (3.17) to 


the corresponding long spatial distance expectation value in hat Minkowski space for a 
massive held. In that case the expectation value decreases exponentially as while in 

de-Sitter we have the oscillatory behavior which is due to the fact that the expansion of 
the universe created pairs of particles. 

For imaginary p, A and A are both real but one of these will give the dominant 
contribution at late times. For the special case of p = 0, we see that both terms in 
(3.17) are similar, but divergent. Taking the limit carefully we hnd a behavior of the 
form W 2 -|- W 2 logzn (see (3.15)). This can be understood as a logarithmic operator [T7] . 
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Namely, we have a pair of operators such that the dilatation operator has the Jordan 

f - l\ 

form D = ( 2 g j. The dilatation operator is not unitary but this is allowed, since the 

three dimensional theory is not arising as the analytic continuation of a Lorentzian 2+1 
dimensional held theory. Of course, for /i > 0, the two eigenvalues are complex. The 
factor rehects the dilution of these massive particles due to the expansion of the universe. 
If we square the wavefunction we will hnd that the probability of hnding the particle goes 
as 1/a^, where a ~ l/r^ is the expansion factor of the universe. More generally, both the 
leading term, as well as the subleading terms in the expansions are related to the quasi 
normal modes for the held 0 in the static patch, see e.g. [HI [19]. The imaginary part 
represents the oscillation of the wavefunction of a massive particle T ~ and the 

real parts represent the dilution factor. The quasi normal modes are | + i/i + n. The n 
comes from the expansion of the wavefunction in power series around r = 0 in the static 
patch, which is in agreement with the dilatation eigenvalues of derivatives of the primary 
operator Oa{x). 


Another important property of (3.17) is the fact that for large /r it is suppressed as 
j-TT/i square this we get which is the Boltzman factor of a massive particle. 


In addition, (3.17) contains a relative phase between the two terms which varies with /i 


between e’’’" for small /i and for large /i. This detailed structure is present for the 
standard Euclidean vacuum. It is also interesting to consider the correlator in the limit 
that the two points are separated by a large timelike distance, with \ x — x '\ = 0 an d 


\r]\ \ri'\. We can get into this regime by replacing w ~ —)■ in (3.17), see (3.15). 

\x\ n 1 r i r 


This extra e*®® gives factors of which imply that the hrst term in (3.17) is now not 
suppressed exponentially for large p, while the second is suppressed by The first 

term is the one expected where we create a particle at t]' and destroy it at r; > r;'. In fact, 
setting 7] = —we see that the hrst term goes as 

It is also interesting to look at the two point function in fourier space. One way to 


obtain it is to Fourier transform (3.14). However, we rederive it by hrst expanding each 


Fourier mode of the held in terms of creation and annihilation operators 


fkiv) 


= fk{v)al + fk{v)a_k 

= , hif,{kri) = {-kr]) = h_i^{kri) 


(3.19) 


fkiv) = i-v)^hi^ikv) , 

h-i^ikv) = -V(e“®®®/cr7) = {hi^ikT])) 




2/i 


where we have also summarized some of the properties of the functions. These creation and 
annihilation operators are the ones appropriate for the early time vacuum, at rj —)• — oo. 
They should not be confused with the out vacuum creation and annihilation operators 
we discussed earlier. The ones in (3.19) are dehned so that the mode functions fk have a 


simple oscillatory behavior at early times, rj ^ —oo. These functions are classical solutions 
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of the equations of motion for each fourier mode. Note that all these equations, (3.19) , 
are true for both real and imaginary /r, here * denotes complex conjugation. As usual, the 
normalization is chosen so that we get the right commutation relations. These can be most 
easily hxed by looking at the early time limit, rj —)■ —oo. The vacuum selection condition 
is given by choosing the solution of the Bessel equation with the right oscillatory behavior 
for early times, or very negative r]. From the small argument expansion of these functions 
we can compute the small rj, r]' behavior of the fourier transform of the correlator 

+ k') [fk{'n)fk{r]')dijl - v) + fk{ri)fk{v'W{v' - v)] (3-20) 


(0fc(h)0-fc(h'))' = 


{vv'V- 

dvr 






k^rjr]' 


+ local (3.21) 


Recall that the prime in (■ ■ ■)' means that we drop the factor of (27r)^(5^(^^ ki). By “local” 
we mean terms that are analytic in k which do not give rise to long distance correlations 
in position space. Notice that the correlator (3.20) depends on the sign of r; — rj'. On the 
other hand, the long distance behavior (3.21), (3.17) does not depend on this sign. This is 
because operators commute outside the lightcone. The terms denoted by “local” in (3.21) 
do depend on the sign of — rj'. Of course we can directly go from (3.17) to (3.21) by 
using the formula 


(Txe 


ik.x 


X 


—2a 




-sUf-o) 


r(a) 


(3.22) 


In order to perform perturbative evaluations of expectation values it is useful to use 
the in-in or Keldysh formalism. In this formalism one performs a time ordered integral 
from the initial time to the time of interest and then one goes back to the initial time in an 
anti-time ordered fashion. Each of the paths corresponds to the perturbative evaluation 
of the bra and the ket respectively. Perturbativative computations reduce to the standard 
prescription but with the path ordering along this special contour. For this purpose we 
need to introduce the correlators with points along different parts of the contour 


(0fc(h)0fc(h'))++ 


fkiv)f_kivyiv - v') + fk{v)fk{vyiv' - v) 

f_k{v)fk{v') 

f_k{v)fk{v') 

fk{v)fk{vyw -v) + fk{v)fkWW{v - v') (3.23) 


Here the -f corresponds to the time ordered part of the contour and the — to the anti-time 
ordered part. 


3.1 Two point functions of fields with spin 

For the case with spin we will not write down the full propagator in de Sitter. Instead we 
will focus only on the late time behavior which is governed by conformal symmetry. For 
an operator with spin it is convenient to contract its indices with a three dimensional null 
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vector e, = 0. This takes care of the fact that we want to look at traceless tensors. It is 
then simple to write the general two point fnnction, in position space, for a general spin s 


(ekO(f)ekO(O)) = 


— 2{e.x){e.x) 


\x 


2A 


X 

with a; = — 
\x\ 


(3.24) 


where e®.0 = Let ns make a few remarks abont the fnnny looking nnmer- 

ator in this expression. This contains the information abont the spin correlation between 
two massive particles. There is a snbgronp of the conformal gronp that preserves the two 
points where the massive particles are sitting. This snbgronp inclndes an SO(3) snbgronp. 
This is easy to see if we view the bonndary as and we pnt the two points at the north 
and sonth pole. Bnt it is still an SO{3) snbronp for any choice of two points. These S'0(3) 
elements are not ordinary rotations in space bnt combinations of varions elements of the 
conformal gronp. Aronnd each of the two points these rotations act like the ordinary ro¬ 
tation gronp. However, the rotation on one of the points is conjngated by reflection along 
the X direction rel ative to the other point. This is related to the term with a factor of 2 in 
the nnmerator of (3.24). Let ns say that x lies along the 3 direction. This implies that if 


one particle has a given eigenvalne nnder the generator J 13 , then the other particle has the 
same eigenvalne nnder J 13 where Jij are the nsnal rotation generators. If one particle has 
some eigenvalne nnder J 12 then the other particle has the opposite eigenvalne. See hgnre 

m 



Figure 4: (a) We see a pair of created particles that approach a given comoving points 

at late time in de-Sitter. There is a subgroup of the conformal group that preserves this 
two points. These include transformations which act as rotations around each of the two 
points as indicated, (b) If two points are separated along the direction 3, then the spins 
of the physical particles in directions, say J 23 are the same. So the two particles would be 
spinning as shown. The spins in the direction J 12 are opposite. 


We now transform (3.24) into Fourier space. We then have a preferred vector which 
is k and we can decompose the operator in terms of the spin around the k axis. More 
precisely, we can say that k lies along the 3 axis, k = (0, 0, k) and choose the polarization 
vectors as 

e = (cosV', sin'^, i) , e = (cosV'', sin'^', —i) (3.25) 
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Note that we only care about e up to an overall scale, therefore we can generically take it to 
be of this form. Then we hnd that the two point function in Fourier space is proportional 
to (see appendix [A|) 


h{e,e, k) 


/2(A,m) 

h{e,e,k) 


oc 



with 


k = 


k 

W\ 


(2s)! 


[s — m) 


(s + m)! 


/2(A,m) 


r(A - 1 + m) r(2 - A + s) 
r(2 - A + m) r(A - 1 + s) 

2 i^i(A - s - 1, -2s, 2 - s - A; 


(3.26) 

(3.27) 

(3.28) 


Here m indicates the angular momentum of the mode around the k direction. We some¬ 
times call this the “helicity” of the mode. The factor in parenthesis in (3.26) is an un¬ 
interesting normalization factor, see appendix (|^. The / 2 (A,m) factor is a phase for 
A = I -|-z/i, with /i real. Note that the hypergeometric function in (3.28) is a polynomial. 

This formula contains an interesting lesson. First, let us recall the formula that gives 
the dimension in terms of the mass for a spin s held 


A± = 



2±V 


(3.29) 


where the last equality dehnes /r. When m = 0 we have a gauge symmetry in the bulk. In 
this case, A+ = 1 -|- s which is the dimension of a conserved current. Let us now consider 
the massive case, but with a small enough mass so that A± are real. The leading late 
time behavior of the expectation value of the held is given in terms of the component with 
dimension A_, which is associated to the more slowly decaying function. An interesting 


feature of the coefficients in (3.28) is that they can change sign. First, let us see this in a 


concrete example. Consider the s = 2 case 




o-2*X 


+ 


4(3-A) . , 6(3-A)(2-A) , 4(3 - A) 


A 


+ 


(A-l)A 


+ 




A 


+ e 


2ix 


(3.30) 

with x = ?/) —-0'. Note that for any integer spin I 2 in (3.27) is a ratio of simple polynomials 
of A. Let us start with A = | where all terms are positive. As A decreases we see that the 
middle term changes sign at A = 1. The fact that it diverges at A = 1 is related to the 
fact that the kinetic term for this mode becomes zero. This is the phenomenon of partial 
masslessness discussed in [22l [23]. For A_ < 1 we have a negative sign. This negative 
sign is a problem because the term corresponds to the expectation value of a held and its 
complex conjugat^ which should be positive. Note that here we are using that for real 


'*’Now that for 0' = 0, e = e* in (3.251. 
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A we can focus on the A_ contribution because it is the largest at late times, forgetting 
about the contribution from A_|_. We should emphasize that the relative sign between the 
various terms is determined by the special conformal symmetry (de Sitter isometries) and 
cannot be changed. In the language of the wavefunction of the universe, we can say that 
the wavefunction is not normalizable. In other words, as positive sign in the two point 
function indicates that the mode has a standard gaussian wavefunction 'h ~ e~‘^ while a 
negative sign implies d' ~ . Note that for the particular case of a massless spin two, 

we get m = 0 and A_ = 0. In that case we see that all components with helicity less than 
two diverge. However, in this case these modes are pure gauge and this divergence is not 
an obvious problem. So the m = 0 case is hne. The conclusion is that very massive spin 
two helds in dS make sense. However, masses within the window 


BAD 


m 

0 < j ^<2 


( 3 . 31 ) 


are not allowed [21]. It remains to be seen if one can make sense of the case m?/H‘^ = 2 
(with A_ = 1). If it made sense, it would be in the context of a theory with an extra gauge 
symmetry, since there is a null state. Recent explorations in this direction can be found in 
|2S|. In this hypothetical theory, the graviton mass would be related to the cosmological 
constant. 


There is one particular example where it is very clear that the range (3.31) should not 
be allowed. Consider a Kaluza Klein reduction with a big internal space. If the internal 
space has a size which is larger than the Hubble radius, then we expect to have massive 


Kaluza Klein gravitons which he in the range (3.31). However, if the internal manifold 


is larger that the Hubble radius, we expect that also the internal dimensions should be 
inflating. 

For general spins the forbidden range is 


m 


BAD ; 0 < ^ < s(s - 1) 


(3.32) 


This can again be found by starting with the case that A_ is real and equal to 3/2. For 
this value all the coefficients in (3.28) have the same sign. As the mass continues to 
decrease, A_ will also decrease, and when A_ crosses one, then one of the terms changes 
sign. Namely, the term with m = 0 associated to the longitudinal polarized held 


/2(A_,m = 0) 


F(A_ - 1)F(2-A_ + s) 
F(2-A_)F(A_ - 1 + s) 


(3.33) 


When A_ crosses one this becomes negative. So, the region A_ < 1 is not allowed, which 
from (3.29) leads to (3.32). As expected, for s = 1 any mass is allowed. 

One might be puzzled by the following question. For large mass we do not expect any 
problem. However, when A = | + f/r, with ji real, the coefficients are complex. If we 
we had a problem with negative coefficients, why don’t we have a problem with complex 
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coefficients?!. In the range where A_ was real, we could argue that the coefficient of 
the dominant power of r] at late times is a hermitian operator whose square should have 
positive expectation value. On the other hand, when we have the oscillating factors 772 =*=*^ 
we cannot argue that the prefactors of each power are hermitian operators. For that reason 
we cannot apply the same argument. The hermitian conjugate of the prefactor of (— 77 ) 2 +*^ 
is the prefactor of (— 77 ) 2 “*/^. But conformal symmetry does not allow a long distance two 
point function between operators with different scaling dimensions. Such a two point 


function would be purely local and it is not governed by (3.28). (3.28) governs the two 


point function of two operators with the same dimension, namely the piece behaving like 
7 ^ 2 +*/^ 77 ' 2 +®^ for example. The analysis in [23] implies that there is no problem in this case, 
see [ 2 S| for a more algebraic discussion. 

Recently the Vasiliev version gravity [271 ESI in de-Sitter has been explored, see e.g. 
[29] . In that case, the higher spin helds are exactly massless and correspond to generalized 
gauge helds. The existence of a forbidden range (3.32) implies that we cannot continuously 


break the higher spin symmetry so that we give a small mass for the higher spin gauge 
helds. This is possible in AdS and it corresponds to dual held theories that are weakly 
coupled. This is saying that the dual to inhation cannot be a weakly coupled large N 
theory (as in [30]). If a dual exists, it should be sufficiently strongly coupled so that the 
anomalous dimensions for the higher spin currents are large enough to pull the bulk higher 
spin helds outside the forbidden range ( 3.32[ ). 

Finally, note that in the limit A —)■ 00 we get l 2 {A,m) = (— 1 )^“™ and 


l 2 {e, e, k) = (e.e — 2e.ke.k)‘ 


(3.34) 


which has the same form as the position space numerator of (3.24), with x ^ k, and the 
discussion around hgure applies. 


4 Three point functions in de Sitter 

In position space, the general form of the correlation function of two scalar operators and 
a spin s operator is given by m 


{0i{xi)02{x2)e\03{x3)) = 


(|x23pd.Ti3 - |Xi3pe.f23)' 


l^isl 


A 3 +S+A 1 —A 2 I ^20 1 ^3H“-5+A2 — Ai 


> 12 ! 


A 1 +A 2 -A+S 


( 4 . 36 ) 


Note that, when s is odd, this function is odd under the exchange 10 2. Therefore, if 
the operators are identical only even spins can arise. 


® In the dS/CFT example for Vasiliev theory [53], the boundary operators defined through the wave- 
function of the universe can sometimes have anomalous dimensions of order 1/7V. However, even in that 


case the bulk higher spin fields are still exactly massless. The argument leading to (3.32) used expectation 
values, not the wavefunction. Also the argument as stated, makes sense only to leading order in G^v 
expansion since we ignored the fluctuations of the background metric. At higher orders we need to be 
more careful about the proper gauge invariant definition of the observable, etc. 
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In the rest of this section we will simply explore these three point fnnctions in Fonrier 
space. For onr eventnal application to cosmology we will be mostly interested in the case 
Ai = A 2 = 3. However, we will also consider the case Ai = A 2 = 2 which corresponds 
to a conformally conpled held in the bulk, or a scalar held with mass = 2H^. We 
will consider this case because it is simpler than the general case but still non-trivial. 
Furthermore, the fonrier transform for the case A = 3 will follow simply from the one with 
A = 2. There have been numerous previous discussions of conformal correlators in Fourier 
space, for a sample see [201 EU EO] 


Given that (4.35) is so simple one can wonder why we complicate our lives to compute 


its Fourier transform. For the cosmological application it is customary to express the cor¬ 
relators in Fourier space, specially since the inhationary universe is also breaking the scale 
invariance. Even from purely theoretical reasons, the work done computing the Fourier 
transforms pays oh when we need to consider conformal blocks, since they essentially 
multiply in Fourier space 


4.1 Two conformally coupled fields and a general scalar 

Here we set Ai = A 2 = 2 and consider a scalar with general A 3 = A. In this case the 
correlator in Fourier space is just a function of fcj = |fcj| which are the lengths of the 
triangle formed by the three momentum vectors. It is a closed triangle due to momentum 
conservation. Scaling implies that it can be a function of only two variables 


Q = 


ki - k2 


k?. 


P = 


ki + k2 

fc,3 


with 


((p(fci)(p(/C2)c^A(fc3))' = ^3 


(4.36) 


(4.37) 


Instead of doing the Fourier transform we will obtain the answer by imposing the 
conformal invariance conditions. When the c orrelator depends only on the magnitude of 
the ki, the special conformal generator ( 2 . 12 ) simplihes and becomes 


b.K = (b.k) 


-2(A 


2)):4 + dl 


(4.38) 


where k = \k\. The invariance condition becomes 

(hi?i + &.i?2 + b.Ki,)k^-^G{q,p) = 0 


(4.39) 


where b.Ki denotes ( 4.38 ) with k —?• ki, which is the momentum of each of the operators. 
For each case we use (4.38) with the corresponding scaling dimension. If we hrst set b so 


that it is orthogonal to k^, and we use the momentum conservation condition, we hnd 

- dlJG{q,p) = 0 ^ dgdpG{q,p) = 0 (4.40) 
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Here we have used that Ai = A 2 = 2. This implies that either G is independent of q or 
p. Since the equation is linear, we can hrst hnd the solution independent of q. Once we 
hnd it, we can hnd the second solution by the analytic continuation ^2 —)■ —^ 2 - Now, we 


assume G is a function of only p, insert it in (4.39), and hnd the equation 


(p^ - l)9pG'(p) + 2p9pG(p) - |A(A - 3) + 2|G(p) = 0 


(4.41) 


This equation has two solutions. We are interested in the one that is regular at p = 1. 
p = 1 is a special kinematic conhguration, it corresponds to a “collapsed” triangle, where 
ks = ki + k 2 , see hgure [^b). Of course, for general triangles we have /ca < ki + k 2 - 
The absence of a singularity for the collapse triangle is related to the adiabatic vacuum 
conditioi]|^ Then the unique solution is proportional to 

G(p) oc 2i^i(2 -A,A-1,1;^^) = 2 Fi{l^ (4.42) 


One interesting property of G{y) is that it is invariant under p —)• —p. This is not the 
case with the whole correlator due to the explicit factor of k^~^ in (4.37). Of course, in 


momentum space, in order to convert a correlator with an operator of dimension A to one 
with an operator of dimension A = 3 — A all we need to do is to multiply it by a suitable 
factor of k. More explicitly 


(0i(fci)02(4)03^(fc3))' = (0i(fci)02(4)03A(4))'fe 


,3-2A 

3 


A = 3 - A 


(4.43) 



Figure 5: (a) Generic conhguration of momenta for the three point function, (b) Nearly 

“collapsed triangle” conhguration with k^ ^ ki + k 2 . 


We now discuss the same result but from the bulk point of view. By a slight abuse 
in notation let us also denote by p and a the bulk helds. Let us further postulate an 

Ht was shown explicitly in ^3^ [3^ [3^ |3^ [3^ that departures from the adiabatic vacuum lead to 
singularities for the collapsed triangle configuration in figure [^b). 
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interaction term in the Lagrangian of the form f Acp^a. The conformally coupled scalar 
has a de-Sitter propagator which is simply obtained from the flat space one 


(‘fk(v)¥’-k(v')y = m 




2k 


(4.44) 


Of course, this can also be obtained from (3.20) setting /i = i/2. 

We now compute {^k^{vo)^k. 2 ^Vo)'^k 3 {Vo))' for a late time rjo. This computation can be 
done using the in-in formalism and it involves two integrals over the interaction lagrangian 


(3pt)'| + 
1+ 


oc 


[fkMh + fk3{vo)I-] 

J —oo V 


i\ 


kik2 


1 

.Vh 


fk3{Vo)G{p) 


(4.45) 

(4.46) 

(4.47) 


We see that the integral depends only on ki -|- k 2 due to the simple form of the propagators 
for the conformally coupled scalars (4.44). The prefactor that depends on fk^ijjo) contains 
both powers of rjQ and gives rise to the two operators with dimensions A± = | ± i/i 


1 




fk^ivo) oc 


(4.48) 


where A and B are some constants. In addition, in (4.48) we see the powers of k^^ ^ 


appearing in (4.37). The factor tIq in (4.47) comes from two factors of t]q one for each 
dimension one operator ip. The factors of in (4.47) are converting the three point 


function involving dimension two operators in square brackets in (4.47), to one appropri¬ 


ate to dimension one operators. Notice that 1 = 3 — 2, so that we can use (4.43). In (4.45) 


we have kept the leading rjo dependence. A more complete expression can be obtained by 


keeping subleading terms from the propagator (4.44), which lead to an additional factor of 
g-«(fci+fc 2 )»?o_ Expanding this exponential gives additional expected contributions from the 


dimensi on tw o part of p. The integral (4.46) can be evaluated as the hypergeometric func¬ 
tion in (4.42). It is interesting to note the early time behavior of hi^{k^ri) oc 2 . 

Then we see that the whole integral behaves as J for early times. As usual 

we make the integral convergent by rotating the contour in to the suitable part of the 
complex ? 7 -plane, r] ~ (positive). An important point is that nothing special happens 
when ks = ki + ^ 2 , therefore the integral is regular there. This is related to the adiabatic 
vacuum condition: there are no particles at short distances. This correlator would be 
singular at this point in an a-vacuum [321 00] for the held a, since we would have both of 
propagator. This also explains why we do not need to consider solutions of 


the conformal invariance equation (4.40) that depend only on q. They have singularities 
at g = 1 and or q = —1. These correspond to collapsed triangles with fci — ^2 T ^3 = 0, 
where the integral representation is perfectly regular. In any case, we also see that if 
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we analytically continue fcs to negative values, then we can indeed find a singularity when 
= fci + /c 2 + ^3 —t 0. This conies from the large \t]\ region of the integral and it represents 
the on shell three point scattering amplitude at high energies. The exponential has some 
T] dependent prefactors that imply that the actual singularity in that region goes as log kt. 
Indeed, we can check that the singularity of G{p) at p = —1 is logarithmic. 

The full expectation value involves also the minus contour which is the complex con¬ 


jugate of (4.46) 


4.2 Three point functions involving two massless scalars and a 
massive scalar 

For our inflationary application we will be particularly interested in considering massless 
scalars, as an approximation to the infiaton. Let us denote the massless scalar by It is 
natural to consider a coupling of the form 




(4.49) 


in the Lagrangian. This is consistent with the shift symmetry, ^ —)■ 4-const ant, which can 
be an approximate symmetry during inflation. 


The mode functions (3.19) for the massless scalar take a simple form 

1 


1 


fkiv) = ^^(1 


ikrj 




—ikrj 


y/2k^ ’ \/2k^ 

Using these mode functions, the integrals producing the three point function are 

/(vov = J 

= [ ^[-klk^T]^-ki.k2{l-ikiT]){l-ik2r])]e"^’'^^’"^^'^a{T],k3] 


= O 


12 


kl2 

Ol2 


= ki + k2 


Tj^ 


- k:i^)Bl„ + + kl- A^|)(l - 

f [{f - p2)(l - p‘)dl + 2(p2 + - 2)(1 - pS,)] 


(4.60) 


(4.51) 


(4.52) 

(4.53) 

(4.54) 


where we used the momentum conservation condition. The main point to note is that the 
results for the massless scalar are obtained by applying a simple differential operator, O 12 , 
to the result for a conformal scalar. It is also important that if we had the integral with 


the other mode functions, the complex conjugates of (4.50), then we would have obtained 


exactly the same operator (4.54). This means that we can pull the operator O 12 out of the 
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in-in contour sum. Notice that now the result also depends on q but in a simple quadratic 
fashion, the complicated dependence is on pQ The q dependence of the Fourier transform 
for general scaling dimensions is more complicated (see the appendix of [22] for an explicit 
expression). Fortunately we will not need it for the inflationary discussion. 

It is also possible to derive this result from the analysis of the conformal invariance 
condition on the three point function. See appendix 

The solution of the conformal invariance condition with the right physical conditions 
is unique, up to an overall constant. On the other hand we could have considered another 
type of bulk interaction such as 

[ eV (4.55) 


instead of (4.49). In this case, an exercise similar to what we did in (4.51) would also give 


us the answer in terms of an operator acting on the result for the conformally coupled 
scalar. However, due to the presence of terms that go like like 1/p^ in the integrand, 
now the operator also involves integrating the result for the conformally coupled scalar. 


Nevertheless, using the equation (4.41) we can check that the result is also proportional 


to (4.54). This can be seen at the level of the lagrangian since (4.55) and (4.49) are the 


same on shell. More manifestly, we can start from (4.49), do some integrations by parts 


and use the equations of motion for the massless and massive scalars to obtain 


j a{VO'‘ = -j V<t«V 5 = -i y V<7V(J^) = t 

where m is the mass of a. 


m 


(VV)e^ = — ai 


(4.56) 


4.3 Two scalars and a massive higher spin state 

We now turn our attention to the case of two scalars with the same mass and a massive 


state with spin. The position space answer is very simple and was written in (4.35). We 


are interested in the fourier transform of (4.35) 


We will focus hrst on the case when the momentum of the operator with spin is much 
smaller than the other two. One naive answer is to set the momentum to zero and write 
an expression in terms of the remaining momenta 


( 0 i(fci) 02 ( 4 )ek 03 (fc 3 ))' . oc \k 


1A1+A2+A-6 






(4.57) 


This would be correct if the Fourier transform was analytic at ^3 = 0. Indeed (4.57) 


captures the analytic part at ^3 = 0. However, the Fourier transform contains also non- 
analytic terms in ^ 3 . These arise from the integration region where the operator with spin 


^When the dimensions are Ai = A 2 = 2 -|- n we get a polynomial of degree n in . 
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is much away than the other two, |x 3 i|, |a; 32 | S> |xi 2 |. Without any loss of generality we 
can set X 2 = 0. Then we can approximate (|4.35[) as 


(0l(fl)02(0)e^03(f3)) 


[xl{e.xi) - 2{e.X3){xi.X3)y 

13^31 I ^11 — A+s 


(4.58) 


k3l>hi| 

We now perform hrst the fourier transform with respect of Xi. This can be expressed as 

1 




oc (V^./ci)^|fc] 


IA 1 +A 2 -A-S -3 


^ |a;i|'^i+^2-A+s 

V = x^e — 2{€.X3 )x3 , where V.V = 0 


(4.59) 

(4.60) 


An important property that we used here is that the vector V is null. This implies that 
all the s factors of ki in {V.kiY are combining into a spin s tensor. We can now perform 
the hnal integral 


\k] 


1 A 1 +A 2 -A-S -3 


3 {V .kiY 


d^X3- 


(4.61) 


|2.g|2A+2s 

This integral has precisely the same form as the integ ral we would need to do to perform 
the Fourier transform of the two point fnnction (3.24), except that e needs to be replaced 
by ki. We can then say that the three point fnnction is given by 


{ 0 Yki) 02 {k 2 )e\ 0 YkY)' - oc \k 


|Ai+A2-A-3| 7 , |2A-3 


fc3-S>0 


1 ^ 3 ! hY-iki^kY 


(4.62) 


where / 2 (F, e,/C 3 ) is the fourier transform of the two point fnnction (3.27 
L — kj 

* \ki\' ^ _ _ 

The fnll small k^ behavior is the sum of (4.57) and (4.62) 
{ 0 Yki) 02 {k 2 )e^. 0 s{h)y . oc \k, 


Recall that 


|Ai+A2-A-3| >, |2A-3 


^3^0 


(4.63) 


where we have not calculated the coefficient in front of each term. 

We will also be interested in some of these correlation fnnctions fnnctions away from 
the squeezed limit. Though the general form seems somewhat complicated, we can focus 
on the component with maximal angular momentum around the k^ direction. In other 
words, we can write the fonrier transform of the three point fnnction as 


{0Yki)02{k2)e^.0Yh))' = [?Yh-k2)r\h 


I s—3 


as{q,p)^ - (4.64) 


where the dots represent terms with lower angular momentum. Th ese ne glected terms 
have the form [e.{ki — k 2 )Y~^{e.kY)'' with z > 0. g, p are dehned as in (4.36). 


® We might worry that in (3.27) e was null, whereas here it is being replaced by the non-null vector 
ki. We can take care of that by projecting onto the spin s component, (fci)®|spin s- 


21 



















When the two scalars have dimensions Ai = A 2 = 2, we find that is a function of 
only p which obeys an equation similar to (4.41). See appendix]^ More precisely, defining 
a^(g,p) = Gs{p) we get 

- l)dlGs{p) + 2(1 + s)pdpGs{p) - [A(A - 3) + 2 - s(s + l)]Gs(p) = 0 (4.65) 

This can be derived using the special conformal invariance condition. 

In the case of massless scalars, when Ai = A 2 = 3, one can also show (see appendix 
(C.2)) that the maximal angular momentum component can be written as 

1 


as{q,p) = Us{p) + ^{p"^ - q^)Gs{p) , 

p / X _ -2sp^Gs + (p^ - 1)(G5 - pdpGs) 

“ A(A - 3) - s{s - 3) 


dpUsip) = -pGsip) 


(4.66) 


where Gs{p) also obeys the equation (4.65). 


5 Four point functions 


5.1 Flat space equal time four point functions 

In order to get some intuition for four point functions, let us first compute some fiat space 
four point functions. We consider four operator insertions at time r/ = 0 and we fourier 
transform in the spatial variables. Let us consider first a massless theory in fiat space. 
The four point function is 

^ ^ 1 A 

(v^fci • • • jr ’ kt = ki + k2 + h + h (5.67) 


Notice that we do not have an energy conservation delta function since we are doing an 
observation at a particular time. We see that we get a singularity when kt —)■ 0. We 
can access this singularity by analytically continuing the momenta. The coefficient of this 
singularity is proportional to the fiat space four point scattering amplitude. A! ~ A. 

We now consider a A<p^ theory. In this case, the contribution of the diagram in figure 
gb), is 




'A 

flat 


where 


A^ 


Akik2k^k4 


J2kitiku,ku,ki) 


±± 


rnat _ 

Ti = 


jflat ^ 2 Ti + T 2 , 


^kt{ki2 + ki){k34 + ki 


rflat rflat rp 

i+- - i_+ - 12 
T 2 = 


2kj{ki2 + ki){k34 + kj 


(5.68) 


45.69) 


ki = ki + k2 , ki2 = ki + k2 , ^34 = ki + k^ , ki = \ki 


(5.70) 
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(a) 


Figure 6 : Diagrams for the computation of the four point function, (a) Diagram for the 
theory, (b) One of the diagrams for the theory. There are two other diagrams that 
arise by permuting the external labels. 


where the ± subindices indicate the sum over all pieces of the in-in contour. Here we see 
a few features. First, we also see a singularity when kt —)■ 0. Its coefficient is again the 
flat space amplitude /[{ki 2 Y — fc|], where we have used that /c^ —)■ 0 to write ^34 = —ki 2 - 


Notice that (5.69) is singular when ki 2 -|- fc/ —)■ 0. This corresponds to setting one of the 


three point vertices on shell, so that energy is conserved at the vertex. Of course, this 
cannot be done for real values of the 3-momenta, but we can reach this point by an analytic 
continuation from kj —)■ —ki, for example. 

The interaction points in hgure should be integrated on both branches of the in- 
in contour, both branches for r] < 0. Alternatively, they can be integrated over all of 
spacetime. For this flat space computation both methods give the same answer, since the 
vacuum expectation value coincides with the vacuum to vacuum transition amplitude in 
the presence of the operator insertions. 


5.2 Four point functions in de Sitter 

Let us hrst consider the case of a conformally coupled scalar. If we have a bulk inter¬ 


action, then the result is essentially the same as (5.67), since the interaction preserves 
the conformal symmetry. The result is 


i^kS^o) ■ ■ ■ ^kSVo)) oc 

i=l 


4 

Vo 


X 


kik2k‘3,ki kt 


(5.71) 


We only needed to multiply by one factor of r/o per external held, where r/o is the time at 
which we compute the correlator in de Sitter spacej^ 


The expression (5.71) gives the leading late time correlators. By putting the operators at different 


times, ? 7 q, z = 1, 2, 3, 4, then one could extract the a subleading term in the t/q expansion which contains 
the part of y with dimension two. 
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As a next step we can consider a A(p^ interaction. Let us first consider the three point 
function 


{^kSvo)^k2(Vo)^k3iVo)y = 


Xril 


8/ciA;2^3 


"Vo 


dr] 


1-00 V 


^iktV _ i 


rvo 


dr] 


1-00 V 


Xr]l 


[ilog{ikt) -i\og{-ikt)] = 


iktv 


'\')o 


-(-jr) (5.72) 


Skik2k3 8fci/c2^3 

The is in the logarithms comes from being careful about the ie prescription for the con¬ 
vergence of the integral for early times. Note also that (5.72) does not have a singularity 
at kt = 0. In this case we have a cancellation between the two terms in (2.10), there is a 
cancellation between the -|- and — part of the in-in contour. Indeed, the hrst hrst term in 
the bracket in (5.72), coming from the 4 - contour, displays the expected singularity. There 
was no such cancellation in (5.71). 


Let us now consider the four point function. This involves integrals of the form 

dr] 


/++ = -/ ^ ^ 


V J- 


VO ,p-iki\v-v'\ 

^lpiki2Vpik34V _ _ 

OO r 27;, 


(5.73) 


where the index indicates the type of in-in contour contribution. For small r]Q, these 
differ from the corresponding flat space integrals only by the factors oil/r] and 1/r]'. These 
factors can be introduced by starting with the flat space result (5.68) and integrating it 
with respect to ki 2 and ^34 




Ti = ^ 


r2 = — 



■ ■ Afc 4 )ds 


7o"A2 , 

Akik2k2,k^ '■ 

1 

Li A ^12 - 

2ki 

r ^ kt 

TT^ 


Wi 



Akik2k^kA 


(=tl)(il)-^±±(^125 ^34) 

Jk34 


+ Li2 


■^34 


-ki iki2 + ki) ,k3^ + ki vr^ 

—> + )- 6 - 


kt = ki2 + k. 


34 


(5.74) 

(5.75) 

(5.76) 

(5.77) 


where we have done the integral. We have only kept the leading r]Q dependence for small 
r]Q which gives us the correlators of the dimension one part of (p. The term 72 , (5.77), is 
a bit subtle. Naively all the contributions from the various pieces of the contour cancel 
out. However, as it happened for the three point function (5.72), we still get a hnite piece 
when we take into account the precise ie prescription for each piece of the contour. 

We see that the result is perfectly regular at ki 2 = kj or ^34 = kj. But we do have 
singularities at ki 2 -|- fc/ —?• 0, ^34 -t- /cj —)■ 0. We also have a singularity at kt —)■ 0. In this 
case 7i, (5.76), goes as 


Ti ~ kt{\ogkt) 




as 


fc —)■ 0 


(5.78) 
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We see that the singularity is very mild because the interaction is relevant and becomes 
small at high energies. Note that the function is still non analytic sX kt = Q because its 
hrst derivative is not hnite at A;* = 0. In general the singularity is produced by 




(6.79) 


where n is the number of particles, which is four in this case, but three in (|5.72|). ^nRa.t 


A' 

Since the 


is the stripped (no energy conservation delta function) flat space amplitudeT 
singularity comes from rj —)■ —cx) we are taking the high energy limit of the amplitude. 



(a) 



(b) 


Figure 7: Kinematic conhguration in the Operator Product Expansion limit or OPE limit, 
(a) In momentum space the momentum ki = ki + k 2 is much smaller than the rest, (b) 
In position space, points form two pairs with the two pairs far from each other. Similar 
limits for higher point functions when a partial sum of the momenta is zero. 




a 


(s) (t) (u) 

Figure 8: Diagrams containing the exchange of a massive held. We consider external 
lines which are massless, or conformally coupled scalars ip. The s diagram is the only 
one that can contribute non-analytic terms in the small kj = ki + k 2 expansion. 


A particularly interesting limit is the A;/ —)■ 0 case. This is interesting because we can 
access it within the space of physical momenta. The non-analytic parts of the correlator 
in this limit determine the long distance correlators in a operator product expansion limit 
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where two points are much closer to each other than they are to the other two pointf^ 
See hgure see [HI [7] . 

When fc/ —)■ 0 we get 




{4pt) 


''tt^ 


8/ciA;2^3^4 


1-2 

kl ^ 12^34 


+ ■■■ 


(5.80) 


The hrst term is what we expect from the operator product expansion with an intermediate 


operator of dimension A = 1. This term comes only from T 2 in (5.75). The second term 


looks like the contribution from an intermediate operator of dimension A = 2. The dots 
denote terms that are analytic in kj at fc/ = 0 or higher order contributions that can be 
interpreted as the contributions of the descendents of the previous operators. Note that 
a term going like kj is n ot an alytic as a function of kj, since it involves a square root. 
However, in the dots in (5.80) there are terms which are independent of kj, which are 


formally bigger than the second term but will not give rise to a contribution in position 
space in the OPE kinematic region of hgure [^b). In this particular example, the second 
term in (5.80) could, in principle, also be a descendent of the hrst but the powers of vr 


suggest o therw ise. Note that 7i does not contribute to the l/k[ term, despite the explicit 
factor in (5.76). Note that the non analytic terms in the small kj expansion can only arise 


from the diagrams in hgure j^s), since there is nothing special happening at small kj in 
the diagrams in hgure [^t),(u). 

The analog computation in Anti-de-Sitter is very similar. But there are two diherences. 
First we only have one of the integrals, similar to /++• Also, the propagator of the inter¬ 
mediate particle is diherent in such a way that the four point function is proportional to 
only Ti- In the OPE limit this produces only a factor of kj which rehects a dimension two 
intermediate state. Indeed a conformally coupled scalar in de-Sitter gives rise to a dimen¬ 
sion two operator on the boundary if we impose the Dirichlet-like boundary conditionj^ 
[ISJIH]. The OPE limit in Anti-de-Sitter was analyzed extensively, starting with |15]. In 
conclusion, in AdS we only see a single operator in this OPE limit, but in de-Sitter we 
see two operators. This will become more clear when we look at the contribution from a 
scalar held with a generic mass. 

We have focused on the singularities in the small kj limit in momentum space. We 
have argued that such terms give rise to interesting power law behavior in position space. 
We are interested in these because they signal the existence of new particles. Let us now 
comment on other contributions to the position space OPE which are less interesting from 
our point of view. In position space there are further contributions to the OPE limit 


Note that this limit looks like an operator product expansion on the three dimensional boundary, but 
it is not an operator product expansion in the four dimensional sense, since the proper distance between 
the operators is always larger than the Hubble radius, even for those that are “near” each other in the 
OPE limit. 

For this mass there is also an alternative choice of boundary condition [32] ■ In that case we get only 
the contribution from the dimension one operator. In other words, only the 1/fc/ term but not the kj 


term in the analog of (5.80). 
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which arise as logarithmic terms, log( 0 :^ 3 /a:^ 2 )) ^^e figure Bb) Such terms arise when we 
do the Fourier transform. If we focus on the regime ki ~ — ^ 2 , then if the large ki limit 
of the four point function goes as l/^i, then we will end up with an integral of the form 
J ^ 2 ) ^ which gives rise to a log |xi 2 | from the large ki integration region. Terms 


that are analyt ic in fcf can indeed give rise to such terms. For example, look at the four 
point function (5.71) which for large ki ~ —^2 contains a factor behaving as 


1 


k\k2kt 


kf ■ 


Indeed we know that the position space version of the four point function (5.71) has the 
form 


z 
z 

(5.81) 

where the variables z and z are defined implicitly by solving the last equations. We see 
that as X 12 — 7 - 0 we get a logarithmic singularity. This can be interpreted as the leading 
correction due to the anomalous dimension of the operator 93 (xi)^ in the theory. This 
is just a simple explicit example, but the fact that such logarithmic terms arise is a generic 
feature of the perturbative AdS or dS diagrams, see |15]. We have seen that these arise from 
the fourier transform from momentum to position space. These terms are less interesting 
for our purposes because they do not signal the existence of new particles. They simply 
reflect interactions between the particles that already exist. Such terms can arise both 
from a contact interaction as in figure ^a) or from analytic terms in kj in the exchange 
diagrams in figure [^t),(u). 


((p(fi) ••-99(^4)) oc 
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5.2.1 Four point function in de Sitter with an intermediate massive scalar 
field 

We can now consider conformally coupled scalars ip coupled to a general massive scalar 
field a with an interaction vertex J Xip'^a. The late time expectation value is given by the 
following expression 

^ Tj^ 2^ 

i^k.ivo) ■ ■ ■ ^kS^o)y = ^ 34 , ki) + other diagrams (5.82) 

7£;(fci2, ^34, fc/) = /++ + /_|— +/—I- + /- 

J —00 '/ J —00 '/ 

where the ± indicate the type of branch in the integration contour. The factor of 2^ 
in the numerator is a symmetry factor. Here {o'kj{v)(^-kj{v')y±± denotes the propagator 
along the appropriate branch, given in ( |3.23[ ). There are two other diagrams and they 
result from the replacements (1, 2, 3,4) —)■ (4,1, 2, 3), (2,3,4,1) respectively. The functions 
I±±{ki 2 , k^i, kj) obey a differential equation that is a consequence of the wave equation for 
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the field a. Namely, from the integral expressions (5.83) one can show that the equation 
(V^ - m^){a{x)a{x'))++ = - x'), or 


- -d, + {kj + 


V 


7]^ 


{(Tkjiv)(^kj{v)y++ = -iv‘^5{ri - T]') 


(5.84) 


implies the equatiorp^ 


1 


kl2 + ^34 
1 


[{ki2 - kj)dl^, + 2ki2dk,2 + - 2 )] /++ = 

[{p'^-l)dl + 2pdp +{m^-2)]G++{p,p') = 

-.1 _ ki 2 , _ ks4^ ^ r< f 

with p= — , p = — , 4± = —G±±(p,p ) 


(5.85) 

(5.86) 


We also have the same equation in terms of the p' variable. We want solutions that are 
symmetric under the exchange p -H- p'. The function I _obey an identical equation to 


(5.85). The functions _and /_^ obey the same equation but with zero in the right hand 


side of (5.85), since (V^ — m'^){a{x)a{x'))+- = 0. 


Using the explicit expressions for the propagator we can either do the integrals (5.83) 


or solve the equations (5.85). The equations contain a bit less information than the explicit 


integral expression. We recover the missing information by imposing the correct boundary 
conditions at the possible singularities of the equation. The possible singularities are 
p = ±1 and p = oo. We demand that the solution is regular at p = 1, and we impose that 
the leading singular behavior at p = — 1 is properly normalized, as we will discuss in more 


detail below, in (5.92). 


For the three point function we had derived the equation (4.41) from the conformal 


invariance condition. This is identical to the homogeneous version of (5.85). For the four 


point function it is also possible to study the constraints coming from conformal invariance. 


If we make the further assumption that the four point function has a structure as in (5.82), 


with a function 4 that depends only ^ 12,44 and kj, then conformal invariance gives us 
the condition 

(Dp - Dp.)G±± = 0 , Dp = ip^- l)dl + 2pdp - 2 (5.87) 


Of course, this implied by (5.85). Is the converse true? It is certainly not, since (5.87) does 


not know about the mass of the intermediate field. We can make a conformal invariant 
function by imposing 


{Dp + m^)G++(p,p') = r(p,p') , {Dp, + m^)G++{p,p) = r{p,p) 


(5.88) 


In deriving this equation we drop some boundary terms at 77 = 0. This is correct for generic masses of 
the intermediate state but might be incorrect for special masses, such as m? = 2 , the conformally coupled 


scalar, or m = 0. Indeed the expressions we will obtain are singular for 77 = i/2 or /r = 3*/2, see (5.931. 
For this special masses it is simpler to do the explicit computation. 
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Here we have included the mass as an arbitrary constant in the left hand side and we 
consider a more general function in the right hand side. For consistency, this function 
should also obey 

{Dp - Dp,)r = 0 (5.89) 

For large m? we can solve for G in terms of a power series expansion of the form 


G++{p,p') = 


r = 


{Dp + rn?) m 


DpV 

r -^ H- 


(5.90) 


This gives a function which is symmetric in p, p', to all orders , thanks to ( |5.89 ). We can 
make r a rational function. The simplest example is to take r = l/{p + p') which is what 


we had above, in (5.85). We could also replace it by r = Dp[l/{p + p')]. This corresponds 
to a bulk theory with an additional contact interaction tuned so as to cancel the term 
coming from a exchange. Proceeding in this way we can add an arbitrary set of contact 
interactions. 

The function G+_ is particularly simple because it obeys the homogeneous version of 
equation (5.85). Then the solution is simply a product of the form 


G-+{P,P') = F{P)F(P') 


(5.91) 


It is this simple product because there is a unique (up to scaling) solution of the homo¬ 
geneous version of the equation (5.85) which is regular at p = 1 and p' = 1. Here the 


function F is the same as the one that appeared in the expression for the three point 


function (4.42). The normalization condition can be set by looking at the region near 
p p' —1. The singular contribution comes from the large rj, p' region of the integral 


in (5.83). In this region, the a propagator simplihes and we can write 


f--\ - —1 








,/2 


2kT 


~ T7“^°s( 1+p)log(l+p') (5.92) 

2k I 


This hxes the function 


G'+-(p,p') = 


TT 


2 cosh (tt/i) 


^ 1 ; + ip, 1 ; 


2 


1 

2 


1 — p'. 


(5.93) 


We have the same expression for G _ 

We were not able to hnd an analytic solution of the equations for the G++ or G _ 

functions for generic masses However, it is easy to derive an approximate expression in 


For nF = 2 which is the conformally coupled scalar, we can check that kjTi in (5.76) obeys the 


equation (5.85). This is an example of an analytic solution for a special mass. Note that 72 is a homoge¬ 
neous solution of the equation for this particular mass, rrp = 2 is a bit special since the derivation of the 


equations (5.85) is not valid due to some boundary terms that come when we integrate by parts to derive 


the equation. For example the function G_|_oc log(p -I- 1) log(p' + 1) does not obey the homogeneous 

equation. This issue does not arise for generic masses. 
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the large p and/or p' limit. For this purpose, let us focus on the one of the integrals, say 
the p integral, and let us compare the integral for the + contour versus the integral with 
the — contour. 


I+± oc [ 




+± > 


I-±(x-j "’'^^'^{a{p)a{p')y_^ (5.94) 


The difference between the two integrals is that we replace ki 2 —)■ —ki 2 in the hrst factor 
and that we consider a different propagator in the second factor. There are situations 
where the propagator will actually be the same. For example, if we consider a regime 
where the two points are spacelike separated in position space, then we expect that the two 
propagators are the same. This means that the non-analytic piece in kj of the propagators 
is the same {akj{p)akj{p')y+±\non-am.iytic = (o-fc 7 (h)o'fcj(h'))-±lnon-anaiytic- See (|3.21[) for the 


explicit expression. This is precisely the region we are interested in in the OPE limit. 
Then we can say 


I+±{ki2, ^ 34 , ki) = -I_±{e^^ki2, ^ 34 , ki) , kj < ki2 


(5.95) 


This equation is true for the non-analytic pieces in kj. Here the precise analytic contin¬ 
uation has been chosen to be consistent with the proper ie prescription for each contour. 
This enables us to compute a part of the /++ or I _integrals from the knowledge of the 


/+_ or /_+ integrals (5.93). 


5.2.2 OPE limit of the four point function with a massive intermediate scalar 

Here we consider the OPE limit, see £gure[^ This involves taking the fc/ —?■ 0 limit of the 
four point function, focusing on the non-analytic terms as a function of kj. In this case, in 
(5.83|) we can replace the correlators {o'^^{p)cr_^^{p'))y^ by their long distance expression 
given in (3.21). This long distance expression is independent of the ± subindices. Inserting 


this into the expressions for the correlator (5.83) we hnd that we can integrate separately 
each of the two terms in (3.21). For each of these two terms, the integrals I±± factorize 
as I±±{ki 2 ,k 34 ,ki) OC J±(/ci 2 ) J±(fc34) witip 


J±(fc42) = r = -(T*)^(fci2)'-^F(A - 1) 

i-oo V 


(5.96) 


[J+(fci 2 ) + J-{ki 2 )] [Mh4) + J-{h4)] = 2(1 + cos(7rA))(fci2fc34) r(A - 1)^(5.97) 

= 2(1 isinh7r/i)(/ci2fc34)“^“*^F(^ -h ipf 


Of course, the expressions in (5.96) are consistent with the analytic continuation (5.95) 


For the second term in (3.21) we get the same as above but with p —>■ —p. Putting all of 


14 


In (3.211 we had factors {rjrj')^ which we should really view as 
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this together, including the numerical factors, we get 




Akik2kzk4 


2'KyJki2kM 


kj 


4 ^ 12^34 


ifi ^ 

(1 + ^sinh7r/i)r(-i/r)^r(- + f/i)^+ 




4 ^ 12^34 


— Ifl 


{1 — i sinh 7 r/i)r(i/i)^r(- — i^Y 


(5.98) 


Note that the region of integration contributing to (5.96) is of order —rj ~ l/fci 2 , ^/k. 


>^34, 


which is much smaller than the characteristic values of rj where the massive propagator 
behaves non-trivially, which is —rj ~ l/kj. We see that in this regime we are very directly 
measuring the long distance propag ator fo r the massive held. Note that ki 2 ~ 2ki ~ 2k2 
as kj —)■ 0. Of course, in expression (5.98) we can replace k 2 ~ fci and k^, ~ /c 4 , /ci 2 ~ 2ki 
and /C 34 ~ 2 /C 3 . 

Notice that (5.98) contains oscillations in the logarithm of the ratio k'j/{ki2k^Y)■ The 


logarithm of this ratio is basically the total number of efolds for which the pair has existed, 
adding the efolds for each member of the pair. Notice that the phase of the oscillations 
is also hxed. It is determined by the mass. Notice that conformal symmetry constrains 


the OPE limit so that we can predict the existence of the two terms in (5.98), but the 


computation of the phase required more information. It required solving the equation for 
the massive held in order to connect the simple boundary condition at early times with the 


late time behavior of the propagator (3.21). The amplitude of the effect goes as e for 


large fi. We can view this leading piece, going as as arising from the interference of 
two processes. One is when we have the usual gaussian evolution for the helds. The other 
is when we create a pair of massive a particles which then decay into (p excitations. Notice 
that the probability to create a pair of massive particles goes as We are getting 

something bigger, of order because we have an interference effect that is sentitive to 
the amplitude and not the square of the amplitude for creating the pair. Therefore, if we 
look at the sky we cannot say for certain where a pair has actually been created, we can 
only detect the deviation in the probability distribution. Notice that indeed the leading 

term in the e~'^^ expansion, comes from the /++ or I _integrals, from the i sinh tt/i factors 

in (5.98). The /+_ or /_ 4 . integrals go like These represent physical processes where 


one produces an actual a particle and it involves the square of the amplitude to produce 
them. 

This effect is very quantum mechanical since we see the oscillations of the wavefunction 
of a quantum particle. Here we are seeing the oscillations because we have the interference 
between producing the pair and not producing the pair of particles. 

We will now repeat the computation in a more indirect way. It is conceptually the 
same as what we did above, but it is easier to generalize to the higher spin case. 
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We now start from the explicit expression for G+_ written in (5.93). The fc/ —?• 0 


limit implies that p, p' —)■ oo. Therefore we simply expand (5.93) for large values of the 


arguments. Each of the hypergeometric functions in (5.93) contains two different powers 
of its large argument 


f (pUi ~ -W + tp) + (2p)''T(v)r(i - tp) 

ZA/TTpa L I I 


+ 


(5.99) 


When we multipy the powers coming from each of the hypergeometric function we obtain 
four terms, only two of such terms give rise to non-analytic dependence in ki. Focusing 
on these terms we hnd 


/+- 


H— 


{4ppy^r{tpyr{^ - ipf + {App'r>^T{-ipfT{^ + ipf (5.100) 


47r(fci2fc34)2 . 

In the small kj limit, we can obtain the results for the other integrals by performing 


the analytic continuation in (5.95). Adding them all up we get again (5.98). 


5.3 Four point function of massless scalars and an intermediate 
massive scalar field 

Here we consider the diagram of the form in hgure |^b) but with massless external helds. 
We denote the massive external helds by We write the coupling as / A(V.^)^cr. As we 
have seen in the section on the three point function, we could get the result for this case 
by applying a differential operator on the result for the four point function of conformally 
coupled scalars. 




A' 


Ap3p3p3p3 

Je = Oi2034^lE{kl2, k34, kj) 


[JE^ki, k 2 , ks, k 4 , kj) + other diagrams] 


(5.101) 

(5.102) 


where /e(A;i 2 , ^ 34 , kj) is exactly the same function we had for the conformally coupled scalar 


in (5.82). The operator O 12 is given in (4.53). Notice that the dependence on (fci — /C 2 ) or 


{ks — /C 4 ) is rather simple and it comes exclusively from the operators O 12 and O34. 

In the small kj limit, the operator O 12 simplihes. ki — k 2 ki 2 = ki + k 2 and the 
operator acts on a power of ku to give 


012^12 ^ -{a - l){a - 2)klf 


Acting in this way on (5.98) we obtain 


JE{kl, k2, /C3, ^4, ki)kj^0 


1 


1287r 


(^12^34)^ X 


(5.103) 
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X 


(1 + isinh7r/i)(- + i/i)2(- + i/i)2r(-i/i)2r(- +i/i)2+ 


l 2 \ ^ 1 

4 ^) ^ ~ ^^^rmr^{- - 


(5.104) 


5.3.1 Four point function with a spinning particle intermediate state 

In this subsection we consider the situation with a spinning particle as the intermediate 
state. We show how the spin information appears in the expression of the four point 
function in the OPE limit. We also derive the precise mass and spin dependence of the 
phase. 

We start from the computation of the /+_ integrals which have the structure of a 
product of two three point functions. More precisely, it is a product of two point functions 
contracted with the inverse two point function of the intermediate massive particle. It 
includes a sum over all the polarization states of the intermediate massive spining particle. 
This has the schematic form 

^+- = y^(Q2Q2QA,m)^ 2 A-.St / a- z{0^_rn0202)' (5.105) 


where the sum runs over all the helicity components of the operator with spin. Here 


J 2 (A,m) is the two point function coefficient for each helicity, see (3.27). The expressions 


for {02020/2,rn) cau be obtained by fourier transforming (4.35) and then projecting onto 


the various helicity components using (3.25). 


Here we are going to focus only on the small kj limit of this expression. More precisely 
we will be interested in the non-analytic terms in the small kj expansion. These are the 
terms that give rise to long distance correlators in position space. For this purpose we 
need the small kj expansion of the three point functions. This expansion involves two 


powers of kj (4.63). For (4.62) we get the same factors of / 2 (A,m) that we got in the two 


point function (3.27). The kj dependence of / 2 (e', e, fc/) comes only through such factors. 
In other words, if were were to set / 2 (A,m) —>■ 1, then we would have no depend ence on 
A)/, all the dependence would be on the angle between e and ki. For this reason (4.57) is 
independent of kj. For small kj we can write the three point function as 


(O2O2OA,m)'fc^_^0 


2^mki2 


1-A 


k]^ ^him) + B,nk 


A-2 

12 


(5.106) 


where Am and Bm depend on the angles only through the angular momentum m projection 
of {kiY on to the kj axis. In other words, they involve the spherical harmonics. Am, Bm oc 
Ysm{d,4>) with cos 6* = ki.kj. At hrst sight one is tempted to keep only the hrst term in 
(5.106) since that is the only non analytic term in kj. But we should be careful, we are 


interested in the non-analytic term of the hnal answer in (5.105), and we do see that we 
have a non-trivial power of fc/ from the inverse two point function in (5.105). Therefore 
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the non-analytic terms in kj in (5.105) are given by 




him) 


(5.107) 


We now see that both terms in (5.106) contribute to the expansion. 



Figure 9: Dehnition of the angles that enter into the answer. 7 is the angle between ki and 
k^. h — h' fhe angle between the projections of ki and k^ onto the plane perpendicular 
to kj. Finally, 6 and 6' are the angles between kj and ki, k^ respectively. 


Now the sum Yhm ^m^-m should be a function of only ki and k^ which is projected 
onto the spin s contribution. Therefore it is the Legendre polynomial of the relative 
angle Ps(cos 7 ). This is not quite what we have in (5.107) since we also have the him) 
contributions. These contributions can be taken into account as follows. Let us hrst dehne 
the angles 


cos 7 = cos 9 cos 6 *' + sin 9 sin 9' cos('0 — h') 

cos 9 = ki-kj , cos9' = k^.ki , COS 7 = ki.k^, 


(5.108) 

(5.109) 


and h ci-nd h' ^ire the angles specifying the directions of ki and k^ in the plane transverse 
to ki. We then expand Ps{cos'j) in dehnite helicity components 


P,(cos 7 ) = (0)P™(P) 


(5.110) 


where the P™ are a version of the associated Legendre polynomials but with a different 
overall normalization. Here we just dehne them as the coefficients of the Fourier expansion 


of the left hand side of (5.110). In other words, we write the left hand side using (5.108) 


and we expand in powers of 
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We now define the functions Q± 

S 

e± = 0±(0,0',V’-^') = XI [^ 2 (A,m)]^' (5.111) 


him) = 


r(| + m + i/i)r(^ + s — ifh) 
r(| + m — ifi)Th + s + ifi) 


(5.112) 


In other words, we dehne these functions by expanding the left hand side of (5.110) in 
powers of ^ and then we multiply each power by the appropriate / 2 (A,m) factor. 

For example, for s = 2 we get 

e, = P,{co.j) - s.„2e.m29 cos(V-V ) +- (3 + i2^)(l + ffl^) - 

We see that for p = 0 we get the usual Legendre polynomial. This is simply due to the 


fact that (5.112) becomes one for p = 0. This is true for all spins. It is also interesting 
to analyze the case of /i —?• oo. In this case we see that (5.112) becomes (—1)®“'". This is 


equivalent to replacing -0' —+ vr, and then changing the overall sign of cos 6' and sin 9' 
we get 

©±1/2^00 ~ .Ps(cosy) , with cosy = — cos^cos^'+ sin6*sin6*'cos('^ —'^') (5.114) 

We see that cosy is the inner product of ki with a reflected version of along the kj 
direction, that flips the sign of only cos —)■ — cos 9'. This is what we expect from the 
discussion around figure 

We also have some helicity independent factors which can be computed by looking at 
the expansion of the highest helicity component. We can write the three point function as 
(see appendix 


{02iki)02ik2)0^^sih))' = [e.ih - k2)Ykf-^-^GM + 


(5.115) 


where the dots are terms of the form [e.{ki — ^ 2 )]* *(/c/.e)* with i > 0. These neglected 
terms do not contribute to the piece with largest angular momentum around the kj axis. 


The function Gg obeys the equation (4.65) . Its solution is 


1 1 1 1 1-p 

Gsip) oc r(- + s + */i)r(- + s-zp) 2 hii- + s + zfi, - + s - zfi,l + s; ) (5.116) 

where p = {ki + ^ 2 )/^/. Here and below the symbol oc means that we have neglected spin 
dependent numerical factors (such as 2, s!, 2®), but we have included the dependence on 
the mass. This is analytic at p = 1 as expected. As p —)■ — 1, it is normalized so that it 
matches the answer expected from the three point vertex 


pV 


*1 


J V 


oc 


(p+ 1)‘ 


(5.117) 
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The large p behavior of (5.116), together with the definition in (5.115) implies that we 
have an expansion as in (5.106) with 


As oc r(- + s + ip)r(-i/i)2 
Bs oc r(i + s - i/i)r(i/i)2*'^ 


(5.118) 


where the ratio Ag/Bg is correctly given by the ratio of the right hand side, with no extra 
factor. 

Putting this together with the angular dependent factors we get 


/ 4 __ oc 


{ki2k 


34. 


kj 


r(-vm- + s + v)^e++ 


-ifj, 


if 


(5.119) 


So far we obtained just one of the integrals. We can obtain the rest of them by using 


the analytic continuation (5.95). We had derived that formula for the case of the scalar 


intermediate state, but we expect the same answer for a spinning particle since it only 


depends on the exponential factor in (|5.94|). This leads to an expression similar to (|5.98|) 

1 


T^\ oc 

4 k34. 


(kuk 


34 ; 


kj 


4 ^ 12^34 


(1 + isinh7rp)r(—i/i)^r(- + s + i/i)^0++ 


k] 


4^:12^34 


-ifj. ^ 

(1 — i sinh7r/i)r(ip,)^r(- + s — z/i)^0_ 


(5.120) 


The angular dependence is all contained within the factors 0± given in (5.111). This I, 


should be inserted into (5.82) in order to get the actual correlator. 


5.4 Four point function of massless scalars and an intermediate 
massive field with spin 

Here we consider external particles which are massless. In this case we have not worked 


out the analog of the operator O 12 in (5.102). However, it is still possible to obtain the 


expression when fc/ —0. The hrst observation is that the angular dependence will be the 
same as in the previous case. The reason is the following. When we consider a three point 
function of two scalar operators and a general operator with spin, then the ratio of the 
sizes of different helicity components is independent of the dimensions of the two scalars, 
in the limit fc/ —>■ 0. This follows easily from the expressions for the fourier transform of 


three point functions that we derived in (4.62) and (4.57). In both of these expressions 


the e dependence does not mix with the dependence on Ai and A 2 . Therefore the only 
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possible difference relative to our discussion for the conformally coupled scalar can only 
come from the overall phase. This can arise due to the different behavior for the maximal 
helicity component. 

We have seen that this maximal helicity component had the form given in (4.66), 


(C.205). For large p the function Gs{p) has the power law behavior 

-..1—A—s I TD ^A—2—s 




+ BsP 


(5.121) 


with As and Bg given in (5.118). Here we used that the normalization condition at p = —1 
is basically the same as for the conformally coupled scalar (5.117), see (C.208). Then the 
power law behavior of the maximal helicity component is given in (C.205) is, for large 
P, 


tts ~ 


Ag A + 1 + s 3 _^_^ Bg A — s — 4 ^_g 


.A “ “ + -r—X- P"""" (5.122) 

4A-3+S 4A-S ^ ^ 

This produces simple extra A dependent factors relative to the answer for the conformally 
coupled scalar. 

Keeping track of these extra factors we can immediatly write the answer for the small 
kj limit 




/C4 / 




I^E\ 


ki<Sik\2 k34 




rkUp 


3 

OC (/Cl2fc34)2X 


(5.123) 


X 


kj 




4^12^34 

kj 

4^12^34 


(1 + isinhvr/i) 


(I + s + ipf 


(I + 5 _ + s + e++ 


— 2/i 


(1 — i sinh TT/i) 


( +g ^/r) ^ ^ 

(f + s + z/r)^ 2 


where the angular dependent factors are the same as before (5.111). For s = 0 this formula 
does not reduce to (5.104). This is simply because here, in (5.123), we have normalized 
the answer to a vertex of the form J which reduces to J for spin 

zero. On the other hand in ( 5.104[ ) we used the verte x f ('V ^)^a. The diference is indeed 
a factor of [(| + i/i)(| — ip)]^ = nv as expected from (4.56) (one factor of per vertex). 


6 Three point function in inflation 

In the previous sections we have considered correlation functions in de-Sitter space. We 
have argued that they are almost completely hxed by the de-Sitter isometries. In this 
section we consider correlators in an inflationary situation. We will assume a standard slow 
roll scenario with an inflaton held 0 which gives rise to curvature huctuations (. In this 
situation, the de-Sitter isometries are slightly broken by the presence of a background time 
derivative 0o- This breaks the scale and special conformal invariance of the correlators. 
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Since 0o is small, this breaking is small and the amount of violation of the symmetries can 
be controlled. In fact, we will see that the symmetries continue to determine the quantities 
that we are interested in. 

In inflation, the fluctuations of the scalar held are well approximated by a massless 
held ((>{t,x) = 4>o{t) + ^{t,x). Then the curvature huctuation ( is given by (see [lO] for 
the precise dehnition of () 



(6.124) 


When a correlator for ^ does not vanish in de-Sitter, then the corresponding correlator for 
C is given by the relation in (6.124). 

The presence of an extra massive scalar held gives rise to a four point function for ^ in 
the de-Sitter case. In the inhationary situation it also gives rise to a three point function. 
We can view this three point function as coming from the same diagram that was giving 
rise to the four point function in de-Sitter, but with one external particles set to (pQ. See 
hgure 10 One important diherence is that now we only have three non-zero momenta to 


play with. In this case kj = and k^ = 0. 



Figure 10: Contribution to an inhationary three point function from an intermediate 
massive particle. The blue line with a cross is an insertion of 0o which is the homogeneous 
time derivative of the inhaton. 


6.1 Intermediate massive scalar field 

Before discussing the general case, we concentrate on the case of a massive intermediate 
held. We can easily analyze the coupling on the right vertex of hgure which comes from 
A J (V0)^(J upon replacing one of the helds 0 —)■ 0o(h) ^md the other 0 —We get 

A f ^[-d^^d^cpoa] —>X f ^[-klr]e^’^^^]d^(j)oO- —t A/cg^o f (6.125) 

Ah Ah Ah 


where 0o = ~v9ri4>o is the proper time derivative. We have assumed that 0o varies slowly 
enough that we can take it out of the integral. The remaining integral has precisely the 
same form as the integrals we had when we were dealing with the coupling to conformally 
coupled helds, except that we set k^ = 0. The vertex on the left side of hgure IT is the 
same as the one we discussed in the de-Sitter case. It can be represented by acting with 
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an operator on the resnlts for the conformally conpled held. In snmmary, the correction 
to the inhationary three point fnnction dne to a massive held has the form 


= MV3)w7j7juOi2lE{ki2, h, h) + Other diagrams (6.126) 

OfC^ fC2 

where lE{ki 2 , ^ 34 , kj) is the fonr point fnnction in the de-Sitter case for external conformally 
conpled scalars, see ( |5.82 ) and O 12 is given in (4.54). Here is the valne of r] snch that 
Vsks ~ 1 , corresponding to the time that the mode k^ crosses the horizon. The other 
diagrams correspond to two other permntations of the 1,2,3 labels. 

We can now obtain the small k^ limit of this expression. This is not identical to the 
small ki limit we took in the previons section becanse now, /C 34 = kj = ks, while before 


we had assnmed ki ^ 34 . Only the diagram explicitly described in (6.126) contribntes 


non-analytic terms in fc/, the others contribnte with analytic terms in the small kj limit. 
As before, it is easiest to write an expression for the /+_ contribntion. This is a special 


case of (5.91), except that now we have p' = 1. In other words, we have 


/+_(p,l) = -F(p)F(l) 
kz 

We can also write a similar (the same) expression for /_4 


limit nsing (5.99). We can also nse 


(6.127) 

We can compnte the large p 


F(l) = 


TT 


-y/2 cosh TT/i 


(6.128) 


For large p, we can hnd I++ and /_by doing the analytic continnation in (^5.95| for the 


p variable. For example, we can get /_(p, 1) = —/+_(e ^'^y, 1). In this way we can hnd 


lE{kl2, ks, /C 3 ) |fc3_5.o — 


1 


TT 


Vksku a/^ cosh Try 

k. 


k?. 


2k 


12 


Ifl 


(1 + i sinh 7 rp)F(—i/i)F(- + ip) + 


2k 


12 


(1 — z sinh 7 rp)F(z/i)F(- — ip) 


(6.129) 


We can now act with the operator O 12 nsing (5.103) to obtain 


07 , 37.3 ) 


9p3p 

rv^ 


ki 




71 


h , 

^ X 


dy^cosh^ vrp \ki 

ks \ (1 — zsinh 7 rp)(| — zp.)(| — ip)T{ip) 

4fci J F(^ + ip) 


X 


(6.130) 
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/ ^3 (1 + ^sinh7r/i)(| + i/i)(| + i/i)r(—i/x) 

\4.ki) r(|-i/i) 


These are the correlators of the scalars, we can obtain the correlators of ^ by nsing 
to obtain 


(6.124) 








ks^O 


c 







2 \2 


eM;,\ 


(6.132) 

(6.133) 


with J defined in (6.131). We have also restored the factors of the Hnbble constant, e is 
the standard slow roll parameter, e = ~ have nsed the standard leading order 

p! 

form for the two point fnnction 


(Cfc C_fc) 


1 


(6.134) 


In (6.132) we have expressed the ratio between the three point fnnction and a prodnct of 
two point fnnctions. It is interesting to compare this to the ratio of two point fnnctions, 
since that ratio appears in the local non-ganssianity which is sometimes nsed a reference 
comparison (in that case the ratio in (6.132) is defined as I/wl)- 


We plot the fnnction J in (6.131) for a couple of valnes of /i in fignre 11 


It shonld be noted that (6.132) represents an additional contribntion to the three 
point fnnction beyond the one present in ordinary inflation HD]. For ordinary single field 
inflation the three point fnnction in this limit is given by the consistency condition and 
gives a constant valne for the ratio in the left hand side of (6.132). However, as explained 


in |13 SH] , this constant valne is simply saying that the effects of the long mode are not 
observable for a short mode observer, any physical effect in the single held theory wonld 
be down by kl/kf [lU]- Notice that the effect we are talking abont is bigger, going like 
[k^/ki]^. It is bigger becanse we now have a second held that leads to physical ehects. 
The fnnction in (6.132) is challenging to measnre becanse of several small factors. First we 
have the slow roll factor e. The fnnction J is also very small. First it is small becanse of the 


dilntion factor Second, as jj, becomes large, there is the exponential snppression 

factor The oscillations are more rapid for larger /r, bnt the overall amplitnde is 

correspondingly smaller. All these snppression factors snggest that one is not going to 
be able to measnre this from CMB measnrements alone, nnless \Mpi is particnlarly large. 
Generically we expect that XMpi ~ 1 or greater. With the 12cm tomography [50] the 
prospects seem to be better, at least in theory. 
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Figure 11: These are plots of the function J(/r, The left side corresponds to /i = 1, 
the right one to /i = 3. The vertical axis displays the numerical value of J. The horizontal 
axis shows log^o(|^) which is the number of decades separating ki and /C 3 . 


6.2 Intermediate massive field with spin 


Here we consider the case of an intermediate massive held with spin. As we have done 
before, we concentrate hrst on the expression for / 4 ._. This expression can be written in 


terms of products of a three and two point function of a form similar to (5.105) 


I -\— — {0^02,0/\^rn=Q)'- 


,2A-3 


^2(0) 


(Oa ^ 3) Inf 


(6.135) 


where only the helicity zero part contributes. This the case because in the two point 
function on the right hand side there is just a single momentum. The three point function 
on the left is the same as the de-Sitter one. The scalar operators have dimension three, 
as appropriate for the massless scalars. The two point function is a new object. It is the 
two point function of a massless scalar ^ and an operator with spin. It would be zero in 
de-Sitter, but it can be non-zero in an inflationary background. This two point function 
is constrained by lorentz symmetry and scaling to be of the form 

{0:,{-k)e\0{k))l, oc ^,k^-\e.ky (6.136) 


Note that only the longitudinal piece of the spinning operator contributes. We have put in 
a factor of (f)Q since we expect it to be proportional to this factor. Note that we are imposing 
scaling symmetry since a constant does preserves the time translation symmetry. Of 
course, at higher orders we will also have a violation of scaling. The equation for violation 
of special conformal symmetry was derived in m 


{h.Ki + h.Ky{Oyh)e^.O{k,))y, = h.dkAOyh)Oykye^.O{ky)', 


/dS 


^ 2=0 


(6.137) 


We can understand this equations as follows. The special conformal symmetry is roughly 
a scaling transformation with a parameter that is linear in T, or the form b.x. A scaling 
transformation is related to a change in y The position dependence of the generator 
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translates into a derivative with respect to momentnm which we see in (6.137). For more 
details see HU¬ 


SO we see that (6.136) is then related to the three point fnnction in a pnre de-Sitter 


backgronnd. The conditions for conformal symmetry imply that the highest helicity term 
in the three point fnnction obeys a simple eqnation. In addition, the varions angnlar 
momentnm components are related to each other. In the end the three point fnnction is 
evalnated at a point where e.{ki — ^ 2 ) = (since ^2 = 0), which picks ont the zero 
helicity component, and g = p = 1, see eqnation (C.186). Since the highest helicity obeys 


the simplest eqnation and the bonndary condition is set at p = —1, we hrst solve the 
eqnation to get the highest helicity component at p = 1 and then we use the constraints 
of conform al syrn metry to go down to the zero helicity component. We are interested in 
looking at (6.137) for the particular case that h.k^ = 0. This implies that only terms of the 
form b.e are relevant in both the left and right hand sides of the equation. This equation 
is analyzed and solved in detail in appendix |C.3[ The hnal result is that the two point 


function can be written in terms of the function Gs{p) that describes the maximal helicity 


part of the three point function (4.66) as 


{OsG.oy = 


r(A-i) 


G.(i) 


T(A - 1 + s) [A(A - 3) - s(s - 3)] 


k^iely 


(6.138) 


and from (5.116) we get 


Gs(l) cx r(^ + s + + s-ip) 


(6.139) 


We now go back to (6.135) and look at the answer for small kj. We only need to 


expand the first factor of (6.135) for large p. The expansion of the hrst factor for large p 


can be obtained by taking the expansion of the maximal angular momentum component 


given in (5.122) and then converting to the zero angular momentum part by using (4.62) 


and (4.57) to obtain 


(03030A,m=o)' OC Ps{cOS 6) 


A 


A + l + s, 
'A - 3 + s 


kt^kl^-^l2(0) + B, 


A 


A 




(6.140) 


with As and given in (5.118). Here cos0 = ki.k^. Together with (6.138) we can write 
the small k^ expansion as 


/+_ OC Ps{cos6){kiky)'i 


3 r(i + s + i/i)^r(| + s — 


■X 
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A(A-3) -s(s-3) 

/c 3 ^ ^ I + ip + s r(^ + ip') r(— ip') 

Wi 


Ua;i 


I — ip — s r(^ + ip + s)T{^ — ip + s) 

+ s r(|-i/i) r(*/i) 

^ + ip — sT{^ — ip + s)T{^ + ip + s) 


(6.141) 


42 

































The result for I _^ is the same. And the /++ and I _contributions can be found in the 

small ^3 limit by using (5.95). Adding these contributions and converting to correlators of 
C we get 


{ah)ah)ah)y 


mki)a-ki)y{ah)a-k3)y 


oc 


ks-i-O 


C = 


-CPs{cos6)Js{iJ-, 

/bl 

0o(h3)-^^ ,^2\2 


(6.142) 

(6.143) 
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r(| + s + i/i)r(| + s — i/i) TT 
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) (1 + i sinh TTjj) 
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+ 71“ (1 — z sinh TT/i) 


I + i/i + s r(—i/i) 

\ — — s r(| — i/i) 

A/x + s r(i/i) 


+ i/i - s r(i + z/x) 


(6.144) 


We have not computed the overall numerical coefficient in ( 6.142[ ). It is simply an s 
dependent numerical constant. It can be computed easily using the techniques in this paper 


once 


chosen. 


the normalization for the higher spin three point interaction A j 


IS 


Note that for s = 0 (6.142) reduces to (6.132) up to a factor of {wrY due to the different 


dehnition of the three point vertex. An important point about (6.142) is the presence of 


the angular dependence through the Legendre polynomial Ps{cos6), where 6 is the angle 
between k^ and ki. See hgure 


12 



(a) 


(b) 


Figure 12: (a) Squeezed momentum contribution where k^ ki, k 2 - Massive helds with 

spin give rise to a dependence on the angle 6 between the vectors k^ and ki ~ —^ 2 . (b) 
In position space, the non analytic terms in k^ give rise to correlators at long distances in 
position space |a; 3 , 2 |, |a; 3 i| S> |a;i 2 |. 
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7 Comments on loop corrections 

A theme in this paper has been the idea that the squeezed limit of three point functions 
or four point functions contains information about the quasinormal mode spectrum of the 
theory in de-Sitter space. Our discussion so far has been limited to tree level effects. In 
this subsection we discuss briefly some aspects of loop corrections. 

One effect of loop corrections is the following shift to the late time behavior of two 
point functions [SQisg (see hgure[I^d)) 

+ ^^i+<57-hM+<5/.) (7.145) 

The term 6fi can be interpreted as a one loop shift in the mass of the held. The qualitatively 
new feature is the presence of a nonzero Sj > 0, which leads to a faster decay at late times. 
This was interpreted in [511 E3] as due to the decay of the particle [M]- In hat space the 
decay might be forbidden by energy conservation, but in de-Sitter it is allowed. The 
particle can decay into two particles of the same mass, and these two particles can fall 
into the cosmological horizon. When all particles have the same mass we get ~ 

[541153] . showing the exponential suppression. Note that the presence of this extra term is 
no problem for the interpretation of each coefficient as an operator in the conformal held 
theory. At three level we had that the two operators were related by A —)■ 3 — A. This 
was a consequence of the fact that the operator obeyed the free held equation in de Sitter. 
Once we include interactions this is no longer the case. We can think of the exponents as 
the quantum mechanical quasinormal mode frequencies, which are complex. They were 
complex at tree level and they continue to be complex now. The only new thing is that 
the correlators are decaying a bit faster than what is expected solely from the dilution due 
to the expansion of the universe. 

7.1 Coupling to a general operator in de-Sitter 

Another interesting loop ehect arises when we have an intermediate loop. For example, 
if the inhaton couples to the Higgs held during inhation, then we can have a coupling of 
the form f (V^)^hh'^. This will then give rise to a loop diagram contribution, see hgure 
[l3|(a). More generally, we can consider a coupling to a general four dimensional operator 
A/ ^2^(4) 'ppig operator could be the pair of Higgs helds h'^h, for example. Then the 
long distance correlator is determined in terms of the three dimensional late time operators 
that the four dimensional operator gives rise to 

0 ( 4 ) ( 7 . 146 ) 

j 

The operators Oaj {x) have dehnite scaling dimensions under three dimensional conformal 
transformations. Note that the operator is an operator in a massive non-conformal 
theory, so it does not need to have a well dehned four dimensional scaling dimension. If 
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Figure 13: (a) Example of a one loop diagram arising from the coupling between the 

inflaton and the Higgs held f (V(/))^h'^h. (b) Example of a diagram where we have a coupling 
to a general four dimensional operator f and we decompose the four dimensional 

operator in terms of three dimensional operators with dehnite three dimensional conformal 
dimensions. The correlator in the OPE regime is governed by these dimensions, (c) When 
we go to an inhationary background we can also compute contributions to the three point 
function, (d) Loop correction to the two point function of a massive held due to a cubic 
self interaction. 


the operator indeed has a dehnite four dimensional scaling dimension, then we can 
say that the smallest Aj is equal to the scaling dimension of the four dimensional operator. 
We then focus on the two point functions of the three dimensional operators 


(OA,(f)OA,(0)) = 


X 


2Ai 


(7.147) 


Using (5.96) (5.97) and (3.22) we obtain 


Ij±± = (±*)(±*) / -ye 


^iki2V 


hE - 


i±± 


IjE - 


kj 


^ 12^34 


-oo ^ 


k^ok: 


f f 

1-00 J 


-v) 




X 


2Ai 


12^34 


87 rsin(7rAj)(l + cos(7rAj))r(Aj — l)^r(2 — 2Aj)cj (7.148) 


Inserting this expression for Ie into (5.82) we get the contribution to the OPE limit from 
the operator Oa^ in (7.146). The full contribution is obtained by snmming over all j 
appearing in (7.146). This is the contribution in the case that the external particles are 
conformaly coupled scalars. 
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The contribution in the case that we have external massless scalars coupled as A / 
is simply given by inserting (7.148) into (5.102) and (5.101). We get a hnal answer of the 
form 




A 2 


Ak\klk] 




k'j 


^ 12^34 


X 


(7.149) 


X 87rsin(7rAj)(l + cos(7rAj))A^(Aj + l)^r(Aj — l)^r(2 — 2Aj)cj (7.150) 


As a more explicit example consider the case that the four dimensional operator is given 
by the square of two massive helds, = cr^, as in the Higgs case. We can read off the 
three dimensional operators that appear by taking the expectation value {cr‘^{ri, x)a^{ri', 0)) 
which is twice the square of the single held expectation value, given in (3.17|) 


2{a{r],x)a{r]\0))^ ~ ^ 


87r® 


rfr] 


\x\ 


3+22/i. 




(7.151) 


+ 


rjT] 


\x\ 


/ \ 3— 


T{^^,rT{--^^,y + 2 


rjT] 


\x\ 




sinh^(27r/i) 

From here we can see that there are three leading operators appearing, with dimensions 

Ai = 3 + 2ifi, Aa = 3 - 2ifi, A 3 = 3 (7.152) 

/ / \ A 

The coefficients Cj are simply the factors of each power of ( ^ j . Here it is clear that 


the coefficient Cj depends on more information than Aj . The dimensions in (7.152) 
arise from the addition of two pairs of particles, the subtraction of two pairs of particles 
or from adding a pair and removing a pair of particles. The coefficients are all suppressed 
by factors of for large /r. 


7.2 Coupling to a general operator in an inflationary background 

We can also consider the case of an inhationary background where the OPE contribution 
arises from a diagram as in hgure [l^c). In this diagram only the operator on the left is 


going to late times and can be expanded as in (7.146). Nevertheless it is still possible to 


hnd a simple expression for the two point correlators of the four dimensional operators in 
this limit 






(7.153) 


(aF — ?7'2 + ie)^^ 

J 

where we have considered a small r] but a hnite rj'. This is determined by de-Sitter 


symmetry plus the assumed late time behavior (7.146). One easy way to understand this 


15t 


For example the third term in (7.1511 has a ^ dependent coefficient while it correspond to a A = 3 


case. 


46 






















is by putting first the operator at some fixed value of x and rj' and then use the de-Sitter 
isometries to move it to a general location The rest of the computation is the same as 


what we did in the case that we were exchanging an intermediate free held. In fact, we do 
not need to do it again, we can simply express the answer we got for a free held in terms 

2 + 


of A = i + iu, and then substitute A —A,-. The hnal result is 




Cfree ( ^ ) 




3. k. 


2 k\ Cfree(‘^_7) 


i^r(^-A)r(A) 


(7.154) 

(7.155) 


with the following dehnitions. C is the same as in (6.133). The Cj were dehned in (7.147) 
and Cfree(A) is the coefficient we get when we had the operator being a single massive 
free held, expressed in terms of the dimension. Finally, the functions J are obtained from 


(6.131) writing 


</(h, ^) = ^(h, 


(7.156) 


l+i^l 


where J is given in (6.131) and J(p, |^) corresponds to the part of J containing the power 


It was noted in [55] that for A* < 2, this is the dominant correction to the three point 
function and that for 2 < Aj there is a bigger contribution of the form kl/k\. However, 
this term has the form of a held redehnition in position space, since it is local in k^. 
Therefore, in principle, we can still unambiguously identify the subleading contribution 
with a 2 < Aj. More physically, we can say that the corrections of the form kl^ jk^^ are 
saying that the statistics of the short modes depend on the features of the long mode at 
the location we are measuring the short mode. On the other hand, a correction which is 
not analytic in kj, implies that the statistics of the short modes depend on the shape of the 
long mode somewhere else, it is not determined uniquely by the local value of the short 
mode. 

Of course, this whole discussion is also valid when the operator is a conformal 
operator. In that case, the leading three dimensional operator also has dimension equal to 
the scaling dimension of the four dimensional operator. Several aspects of this situation 
were discussed in [55]. Other simple examples include the Yang Mills held strength, with 
A = 4 or pairs of Standard Model fermions with A = 3 if they are massless. 


^®More explicitly, set 7=0, then (7.146) plus the dilatation symmetry implies that the answer is 
We can then use a special conformal transformation which moves the bulk point but leaves the 
position of the boundary point fixed to get the general answer. 
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8 Conclusions 


Inflationary cosmology provides us with a natural high energy accelerator. The late uni¬ 
verse represents the detector output of this accelerator. The presence of new particles, 
beyond the inflaton, lead to subtle imprints on the cosmological primordial fluctuations. 
Of course, the main features in these fluctuations are set by the inflaton. The fact that we 
can approximate inflaton fluctuations by a nearly free held implies that these huctuations 
are very close to gaussian. Self interactions of the inhaton, or interactions between the 
inhaton and other particles give rise to non-gaussianities. In fourier space, these non- 
gaussianities are a function of the momenta of the huctuations. When we look at the 
simplest three point function, we have a function of a triangle. When this triangle has 
a very small side it is interesting to look at the behavior of the correlator as a function 
of the ratio between the small side and the other two large sizes k^/ki. As a function of 
this ratio the non-gaussianity displays a characteristic power law behavior which signals 
the presence of new particles H El m E]. This power law behavior is a generic feature of 
correlation functions for quantum held theory in de-Sitter space. The particular powers 
that appear seem to be connected to the quasinormal modes for the quantum held theory 
in the static coordinate patch. We have not given a proof of this fact, but it is true in the 
examples we considered and it seems to be a reasonable expectation. 

Massive particles with masses of order H give rise to interesting signatures in the non- 
gaussianities. These signatures arise for very reasonable couplings which we generically 
expect to be present. They contain information about the mass and the spin of the particle. 
When m > ^ there are oscillations as a function of the ratio between the scales. These 
oscillations arise due to the quantum mechanical nature of the problem and the fact that 
we are doing an interference experiment. It is the “cosmological double slit experiment”. 
The oscillation frequency tells us about the mass. The phase of the oscillation is also hxed 
in terms of the mass |E] • Other oscillatory effects were discussed in [56l |57] and they differ 
because the oscillations are due to oscillatory features in the potential which are in phase 
with the background evolution. In other words, in those models the oscillations depend 
on the value of k while here they depend only on the ratio between two momenta. In 
principle, with enough data, one can distinguish between these two types of oscillatory 
behavior. 

In terms of measurability, if the couplings are Planck suppressed, then it seems impos¬ 
sible to measure this through the CMB or large scale structure. (See e.g. [5H1 EHl EHl El] for 
a discussion of measuring these effects via large scale structure.) But it might be possible 
using the 21cm tomography [50] . It would be nice to analyze what range of masses and 
couplings are detectable. 

We should stress that the effects we have discussed in this paper are not the largest non- 
gaussianities. They are the effects that most clearly reveal the presence of new particles. 
It is likely that the non-gaussianity will be discovered hrst, and only then, a more precise 
study will reveal the new particles. 

In this paper we have worked in the approximation of a quantum held theory in de- 
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Sitter space. In this case the conformal symmetry is exact and its action is very clear. If 
we consider a theory with dynamical gravity, the situation is slightly less clear because 
local operators are not proper gauge invariant observables. Also the correlators are not 
de-Sitter invariant in the sense we discussed here [51]. This is due to the fact that the 
graviton propagator cannot be chosen to be de-Sitter invariant. This lack of invariance 
should cancel out in proper gauge invariant observables. But the standard four point 
function we have been considering is not a proper gauge invariant observable. It was 
shown in [511 how to modify the action for the conformal generator to make the four point 


function conformal invariant. It would be nice to see if a statement like (1.1) can make 
sense to higher orders in the gravitational perturbation theory. 

It would also be nice to see how new particles affect correlation functions with external 
gravitons. In other words, we could also use non gaussianities that involve tensor modes 
in order to detect new particles. 
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A Derivation of the two point function in momentum 
space 

Let us write down the special conformal generator in momentum space. Its action on a 


scalar operator was written in (2.12). To describe its action on an operator with spin. 


it is conventient to introduce a null vector T, P.P = 0 and to contract the operator with 
this vector: e^.Og = Here is a traceless symmetric tensor. Our 

choice of e produces also a traceless symmetric tensor when we take s factors of e. Then 
the action of the special conformal generator is given by 

b.K = -(A - 3)2b.dk - [(b.k){dl) - 2{k.dk)(b.dk)] + 2[{€.dk)(b.d,) - {b.e){dkA)] (A.157) 

When we take the derivatives with respect to e we can treat its three components as 
independent. The error we make by not taking into acount the = 0 condition cancels 
out at the end. This is easily seen since the last bracket vanishes when acting on a term 
of the form e^. 

We now consider a two point function and we make an ansatz of the form 



2=0 


(A.158) 

s 

III 

« = A 

k • 

p = te 

(A.159) 


Note that h is dimensionless. When we act with the conformal generator we think of 


the momentum in (A.158) as that of the hrst operator. When we act with the conformal 
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generator (A.157) we can move through the factor of so that now the operator only 

acts on h. Note that the part of the operator acting purely on vanishes due to the 

usual conformal invariance condition for scalar operators. We only have to worry about 
the case where only one of the derivatives of the operator (A. 157) acts on k^^~^ and the 


other acts on h, as well as the case where both act on h. Note that the term k.dk measures 
the total power of k and it can be commuted past the h.dk, producing a —{h.dk) term in 
the process. We then use that k.dk produces a simple factor of (2A — 3). 

We will also assume that we choose b in such a way that 


b.k = 0 , b.e = 0 , he 7^ 0 


(A.160) 


This condition implies that the term with (b.k) in (A.157) does not contribute. We then 


obtain the following expression for the conformal invariance equation 


(A - l){b.dk) + (2A - 3)^(ha:) + {e.dk){ld,) 


h = 0 


(A.161) 


Here we have used (. 

4.160 

) in order to drop various terms, includding the last term in the 

spin factor in ( 

A.157 

’). Now we use the following properties of a,l3,p. 


dj:a = 7(7*— ka) , 

def3 = 0 , 


k 

= k , 


dk/3 =-{e - k(3) , = 0 


dep = e 


Using these properties it is possible to show that (A. 157) becomes 


0 = [(A - l)di3 + (2A - 3)adp + + p)9^]c?p] h 

Writing h = {a/3)^F{-^) we hud that the above equation becomes 

{y - 2)yF" + [(A - s)(2 - y) - 1]F' + s(A - 1)F = 0 ^ 


F= 2Fi(A-1,-s,A-s-^;|) , 


y = 




(A.162) 

(A.163) 

(A.164) 

(A.165) 

(A.166) 


This equation is of the hypergometric form and we are interested in the solution that is 
a polynomial for integer Writing e and e as in (3.25), with k = (0, 0, A;), we hnd that 
y = 1 + cos('0 — 'll)'). Inserting this in (A.166) we can see that it is proportional to 






(-l)"-™(2s)! r(A-l-m)r(A-l + m) 
(s — m)!(s + m)\ r(A — s — 1) r(A + s — 1) 


= 2 i^i(A-s-l,- 2 s, 2 -s-A;-e*^) 


(A.167) 

(A.168) 


It is the Jacobi polynomial ® — 
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(A.169) 


with 

In position space the two point fnnction has a very simple expression 


{2pt) 


[X^.X2{VvV2) - (X,.V2)(X2 .Vi)Y 
(Xi.X2)^+* 

[0:12(61.62) - 2(ei.a:i2)(62.a:i2)]^ 


|a^l2| 


2A+2s 


(A.170) 


A = (A+,A-,A) = (1,-f^f) , 0 = X+X- + X\ 1/= (0,-2e.x, e) 
where the hrst expresion is written in embedding coordinates, snmmarized in the last line, 


with V.X = = 0. We have checked in some cases that the Fonrier transform of (A.170) 


does indeed match (A. 166). 


Let us record here the momentum space expressions for s = 1 


3-*x + 2^^^ + 


A - 1 


(A.171) 


The expression for s = 2 was given in (3.30) 


When A = I we see that the phases simplify and we get a result that the angular 
dependence goes as h oc We think that this is due to the fact that in this case, we 

are getting that the prefactor —)■ 1 so that the Fourier transform is completely local. 

In this case we expect that the non-local parts in x should arise from a term that has a 
log k in momentum space. The simple angular dependence arises because the local term 
should have a simple angular dependence. 

Another interesting limit is A —)■ cx). This corresponds to the massive case. Note that 


since (A. 168) is a rational function of A, it does not matter in which direction in the 


complex plane we take A to inhnity. Then the result becomes 

h oc (F.e — 2 e.ke.ky oc (cos('0 — y') — 1)^ 


(A.172) 


The combination that appears here is similar to the one appearing in the numerator in 
position space, 

(F.e — 2e.xe.x)^ (A.173) 


see (|3.24|), where x = This combination is related to the following fact. If we have 


two points on the boundary. Then they are left invariant by an 50(3) subgroup of the 
conformal group. That 50(3) subgroup acts like ordinary rotations around each of the 
two points. However, one of the rotations is conjugated by a reflection along the line that 
joins the two points. Notice that 


e.e — 2 e.x e.x = e.Rxe 

where Rx is a reflection along the x direction. 


(A.174) 
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B Three point function with a pair of A = 3 operators 
and a general scalar operator 

In this appendix we obtain the expression for the correlator of two scalar operators with 
dimension Ai = A 2 = 3 and one scalar operator with general dimension A by using the 
special conformal invariance conditions in momentum space. Using (4.38) we can write 

0 = + (6i3)[-2(A - 2)-!-4, + I A (B.175) 


i=l 


Picking first b so that b.k^ = 0 and using the momentum conservation condition, we end 
up with the equation 


1 - 2 ^ 4 . + 2 ^ 4 , + (dl - dl)]F = 0 


fel - ^k21 

2 


0 = [2{qdp - pdq) + {p^ - q^)dqdp]F = 0 , 


q = 


h - k2 
k?, 


p = 


ki + k2 

k?, 


Defining Q = q^ and P = p^, then we find that (B.177) becomes 

0 = [{dp - Oq) + {P- Q)dpdQ]F = 0 
We take then dq of the above and get 

{-l + {P-Q)dp)dlF = 0 


(B.176) 

(B.177) 


(B.178) 


(B.179) 


Then this implies that F = Uo{P) + Qui{P) + Pwi{Q) + Wq{Q). In other words, F is 
a polynomial of up to degree 1 in Q or Inserting this now into ( |B.178[ ) we get that 
the variables separate and we get 


dpUQ + {Pdp - 1)mi = a = dqWQ + {Qdq - l)wi 


(B.180) 


We will be interested in solutions with a simple Q dependence. Thus we set Wi = Wq = 
0 = A. Then we see that we can determine Uq from Ui. We can also write 


F = U{p) + ^{p^ -q^)G{p) , 


U'{p) = -pG{p) 


(B.181) 


We can now take b along ks in (B.175) to obtain 


0 = -{A - 3)A{p‘^ - q‘^)F + 2q{-2 - p'^ + q‘^)dqF + {p"^ - q‘^){q‘^ - l)d‘^F+ 


If we had assumed that the operators at 1 and 2 both have dimensions A^, we would have obtained 
that it is a polynomial of degree Ag — 2 in Q, when this quantity is an integer. Otherwise, we get a more 
complicated function. 
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+2p(2 + q^ - p^)dpF + 2pq{p^ - q^)dqdpF + {p^ - q^){p^ - l)5pF (B.182) 


We can now insert the expansion ( |B.181 ). We then find a polynomial in q whose highest 
coefficient is once we nse the eqnation for the derivative of U in (B.181). This highest 
coefficient gives the eqnation 


- l)G'\p) + 2pG'{p) - (A^ - 3A + 2)G{p) = 0 


(B.183) 


which is the same as the eqnation for the Ai = A 2 = 2 case (4.41). The snbleading 
coefficient gives 

U(p) = (B.184) 


solve (B.183) and determine U from (B.184). 


This eqnation is compatible with (B.181) and (B.183). Thns this the final solntion, we 


In the particnlar case of the (VC)^o' interaction, we have seen that we can get the 


answer by acting with the operator in (4.54) on the resnlt for a conformally conpled scalar 


with an J ip'^a interaction. This then prodnces the expression 

G{p) = -2/ + 2pf + - l)f = A(A - 3)/ , U{p) = (p2 - 1 )(/ - pf) (B.185) 

where / = Gconf /2 where Gconf is the Fourier transform when the two scalar operators have 
A = 2 (see ( 4.41[ )). This is compatible with (B.184) and it gives the properly normalized 
interaction with the vertex in (4.56). 


C Thee point functions of two scalars and a general 
spin operator 

In this section we present the derivation of the three point functions with two scalar 
operators with special dimensions and a general operator with spin with general dimension. 
Of course, the expressions of the three point functions are simple in position space, see 


(4.35). We present the derivation in momentum space, which, in particular is giving us 


the expressions for the Fourier transform of (4.35) 
We consider the correlator 


(Oa^Oa^GOa) = + 


(C.186) 


i=0 


q = 


\ki\ - 1 ^ 2 ! 
1^1 + ^2! 


7 = e-(^i-4), 5 = e.{ki + k2) 

\ki + K 2 I 


Since we have written all momenta in terms of fciand k 2 , we can forget about the action 
of the special conformal generator on the third operator. We now consider a vector b that 
obeys 

b.ki = 0 , b.k2 = 0 , b.e 7 ^ 0 (C.187) 
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Then the action of the conformal generator can be written as 

b.Ki + h.K2 = (b.dki + bJk^){-2{Ae - 2^+ ki.dk^ + h^dk^) + (b.di^^ - b.dk^){ki.dk^ - k2.dk^) 
= (A - 2)(6.4i + (C.188) 


where we used the overall scaling dimension of (C.186). We will need the following equa¬ 
tions 


{b_.dk^ - = 2{e.b) _^{b.dk^ - b.dk^)6 =_0 

(b-dki + b-dk^h = 0 > (b.dk^ + b.dk^JS^ 2{e.b) 

{ki.dk^ - k2.dk^)'^ = S , (ki.dk, - k2.dk2)b = 7 
(h.dk, - k2.dk,,)\ki + k2\ = \ki + k2\qp^ 

(ki.dk, - k 2 .dk^)q = p{l - q^) , (h-dk, - k 2 .dk^)p = q(l - p^) (C.189) 


We then obtain 


0 = {d.y[p(l-q‘^)dq + q(l-p'^)dp + (A + 2(A,,- 3) - s)qp + 6d.y +'yds] + (A-2)ds} X 


X 


J2lW-^^.(q,p) 


(C.190) 


1=0 


which implies the following recursion relation for the a, coefficients 

ai-i = - ^ -2 + i) ~ ^ ~ + (^ + - 3) - s)qp} a^+ 

+ (i + l)ai+i] (C.191) 

This recursion relation determines all the a* coefficients once we know as{q,p). (We should 


set Os+i = 0 in (C.191)). Notice that corresponds to the highest helicity component 


around the k^ axis. The crucial point is that conformal symmetry relates the various 
helicity components. 

Note that we have not yet imposed all the special conformal invariance conditions, 


only the subset in (|C.187|). The other conditions imply equations that determine Ug. 

(C.192) 


These equations can be most easily obtained by picking b so that it obeys 

b.e = 0 


This still leaves two independent components for b. In o rder to write the remaining equa¬ 
tions is it convenient to go back to the ansatz (C.186) and replace ki + k 2 = —k^. In 
that way, the coefficients depend only the \ki\, i = 1,2,3 and the only dependence on 
the vectors is in 7 and 5. One then hnds that the e-independent part of the conformal 


generators (A. 157) does not act on 7 or h (there is a cancellation between different terms 


in (A.157)). Finally, the e dependent piece in (A.157) has a simple action. The result is 
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that we get two equations, one per b component that obeys (C.192), of the form 


0 = 


0 = 


hKi-K 2 )- 26 d^^dk, 
z fcs 


^j^^^A+2iA.-3)-sJ2YS^-^a,{q,p) (C.193) 


i=0 


Ks--{Ki + K2) + 26ds-dk, 
z fca 

ki = -2(A. - 2)-t4. + dl , 

rCj 


\h\A+"^A.-3)-sJ2Y5^-^a,{q,p) (C.194) 


i=0 


Q = 


h - h 

k?. 


ki + k2 , , 

P = —- C.195 

kp. 


Here Ki is the same operator we encountered when we considered the action of special 
conformal generators on functions that depended only on ki = \ki\. The 9^. derivatives in 


(C.193), (C.194), by dehnition, do not act on 7 , 5. An important property of the equations 


(C.193), (C.194) is that they do not mix the highest helicity component as with the rest 


of the functions. Therefore we get two closed equations for the highest helicity component 
which are similar to the equations we would get for a scalar correlator, but different in 
detail because of the extra power of in the above ansatz. The two equations are 


0 

0 


[ 2 (Ae - 2){qdp - pdq) + - q^)dqdp] as{q,p) (C.196) 

{-[(A - 3)A - AAl + 18Ae + 4AeS - s{s + 9) - 18](Ae - 3){p^ - q^)- 
2q'[2Ae — 4 — (1 + s')p^ T (5 — 2Ae + s)(f‘\dq-\- 
+ 2 p[ 2 Ae - 4 - (1 + s)q^ + (5 - 2 Ae + s)p^]dp 


+(p 2 - - l)dl + - q^){p^ - l)dl] a,(g,p) 


(C.197) 


We will now discuss these equations in turn for two special cases, Ag = 2, 3, which 
corresponds to the conformally coupled and massless cases respectively. 


It will be interesting for us to solve the recursion relations (C.191) at g = p = 1. 


This is a point where the aj(g,p) functions are required to be analytic. Then we see that 
the terms with derivatives drop out since they have a factor of (1 — q^) or (1 — p^). We 


can simply solve the recursion relations by going back to (C.190), setting q = p = 1 and 


then making an ansatz of the form * 01 ( 1 ,1) = Replacing this in the 

equation (C.190) we hnd 


S 


^ 1) = 5'*(1 - 2 Ri(A + Ae - 3, -s + Ae - 2 , A^ - 1; (C.198) 

i=0 


C.l Two = 2 operators and a general operator with spin 


In this subsection we set Ag = 2 in the above discussion. The equation (C.196) becomes 


very simple, and it says that depends only on p or only on g, as we will see below we 
want the solution that depends only on p, as{p)- Inserting this into the second equation 


( IC.197D , we hnd 

0 = [(p^ - 1)9^ + 2(1 + s)pdp - (A - 1 + s)(A 


0 ] (^s{p) 


(C.199) 
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which reduces to (4.41) when s = 0. The solution that depends only on q is set to zero by 


demanding regularity at g = ±1. For large p this solutions behave as 


as{p) ~ p 


l-A-s 


.,A-2-s 


p 


(C.200) 


Inserting this solution into the recursion relations (C.191) we get all the other helicity 


components. For example, for massive spin two we get 

2q 

[(/ - l)dpa2{p) - Apa2 


(C.201) 


ao{q,p) = 


[(1 + (3 + 2A)g^)p(p^ — l)dpa2+ 


A(A-l) 

+{A(p2 - 1 ) + + 4A - A^ - 4p2 - 5Ap2 + 2AV)}a2] (C.202) 

Here we have used the equation for 02 , ( |C.199 ), to simplify the expressions. In general, 
we find that is a polynomial of degree s — i in q. Note that he singularity at A = 1 is 


the same as the one we discussed around (3.30) 
For the s = 1 case we get 


Uq — 


g[(l - p^)dp + (A - 3)p]ai(p) 
A - 1 


(C.203) 


Note that for odd spin, the full answer is odd under the 1 -f-)- 2 exchange. So that if the 
operators 0{ki) and 0 {k 2 ) are identical the three point function vanishes identically. 

Let us go back to general s. For g = p = 1 we can solve the recursion relations as in 


(C.198). Specializing also to the case that 7 = <5 we get that 




i=0 


F(A-l + s) 


(C.204) 


This relation is important if we want to to set ^2 = 0 and ki = k^. This sets j = 6 and in 


addition. Then the relation (C.204) is specifying the value of the zero angular momentum 


component in terms of the maximal solution as{p) evaluated at one. This relation is useful 


because as{p) obeys the simple equation (C.199) 


It is possible to solve the recursion relations explicitly for large p. For large p, has 


two possible power law behaviors, coming from the large p solution of (C.199), see (5.121). 


Inserting each of these two possible powers in the recursion relations (C.191) we obtain 


results consistent with (4.63), where each term comes from each of the two possible powers 
at large p. 


C.2 Two Ag = 3 scalars and a general operator with spin 


Here we set Ag = 3 in the discussion of appen dix [C| The first equation (C.196) becomes 
identical to (B.177[) whose solution is given in (B.181) 


a,(g,p) = Us{p) + - q^)Gs{p) , 


dpUsip) = -pGs{p) 


(C.205) 
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Inserting (C.205) into (C.197) we get a couple of equations. The first is that Gs{p) obeys 
the same equation as in (C.199) The second equation gives Ug in terms Gg 


Usip) = 


-2sp^Gg + - 1)((T^ - pdpGg 

A(A-3) -s(s-3) 


(C.206) 


which for s = 0 reduces to (B.184). Again, inserting these solutions in the recursion 


relations we get the rest of the ai{q,p). 

We now consider the solution of these recursion relations for q = p = 1. This is given 


by (C.198). Setting Ae = 3 and 7 = h in (C.198) we hnd that it vanishes X]i=o 1) ~ 0. 
However, for our application, we need to take the derivative of the expression ( |C.198 ) with 
respect to — ds before setting them equal. This gives a nonzero answer. We hnally get 




2=0 


7=(5=1 


s! 2 T(A) , , 

s!27A-l)r(A-l) 


(C.207) 


-2s)a(l) 


r(A-l + s) [A(A - 3) - s(s - 3)] 


where we have also used (C.205) and (C.206) at g = p = 1 (using that Gg{p) is analytic 
at p = 1 ). 

It is also interesting to consider the case p = — 1 which corresponds to the singularity 
at ki 2 = —kj. In this case we expect the singularity to be of the form 


-t kik 2 J —e 




OC 


{P + 1 )^ 


(C.208) 


where we focused on the most singular piece. We see from (C.205) that this will come 


from the second term in Ug, since the U term involves an integral that will reduce the order 
t p : 

- 1 . 


of the singularity at p = —1. This imposes the same condition on Gg as (|5.117|), namely 
(p+Tp ^ 


G 


C.3 Calculating the two point function in inflation 

In this subsection we calculate the two point function of a dimension three scalar operator 
(associated to a massless held) and a general dimension A operator with spin. We begin 
with the identity 

(6.i?)(03(-fc)ek0(fc))U OC bA,{03{h)0s{k2)G0{k)y\-,^^, (C.209) 

where k = —(fci + ^ 2 ). We consider b so that b.k = 0. This will be enough to determine 


the two point function. The ansatz for the two point function was given in (6.136). The 
left hand side of ( |C.209 ) is 


{b.K){Os{-k)G.O{k))[^^ = -2(A - l)k^-^{e.kr-\e.b) 


(C.210) 
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We have only one generator acting on the left hand side since we can write all the explicit 
momenta of the two point function in terms of the momentum of the second operator, by 
using the momentum conservation equation. We can now compute the right hand side of 


(C.209). The hrst observation is that if the derivative acts on a scalar, such as fc,, then 


by the time we set ^2 —>■ 0 the direction will become proportional to k, which contracted 
with b will give zero due to our condition on b. Then the only terms will come from the 
derivatives of terms that contain 6 .^ 2 - This will give rise to a derivative of the form d^ — ds- 


This gives us a result proportional to (C.207). Equating this to (C.210) gives us 




r(A-l) 


G.(l) 


r(A - 1 + s) [A(A - 3) - s(s - 3)] 


k^{e.ky 


(C.211) 
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